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ABSTRACT 

The  analysis  of  the  solidification  curve  for  the  traditional  cases 
of  equil  Ibriuin  and  normal  non-equilibrium  solidification  has  been  extended 
to  include  computation  for  a  binary  alloy  of  any  phase  diagram.  An  addi¬ 
tional  case,  taking  into  account  diffusion  within  the  solid  phase,  has 
been  analyzed  to  give  a  solidification  curve  sensitive  to  the  mode  of 
solidification,  in  particular  the  size  of  the  solidification  element. 
Expressions  have  also  been  derived  for  each  of  the  three  solidification 
cases  in  order  to  compute  the  rate  of  solid  formation,  the  interface  and 
average  compositions  of  each  phase,  the  rate  of  change  of  the  composition 
with  the  progress  of  the  solidification,  and  the  relation  of  the  amount 
of  non-equilibrium  phase  to  the  size  of  the  solidification  element. 

Numerical  analysis  was  used  to  calculate  the  change  of  the  concentration 
profile  on  a  microscale  within  a  dendrite,  due  to  diffusion  within  the 
solid. 

Apparatus  and  control  equipment  were  designed  and  assembled  for 
controlled  Ingot  solidification.  In  addition  to  the  extended  solidifica¬ 
tion  heats  described  below,  the  equipunent  permits  (a)  interrupted 
solidification,  (b)  isothermal  solidification,  (c)  programmed  solidifica¬ 
tion,  (d)  unidirectional  solidification,  (e)  solidification  under  vibration, 
and  (f)  solidification  with  stirring. 


ii. 


Ingots  of  •luminu*~4, 5  p«r  cent  copper  were  solidified  over  periods 
ranging  front  twelve  to  one  thousand  hours.  The  resultt  demonstrated  that 
the  increase  in  the  size  of  the  solidification  element  concomitant  with  a 
decrease  in  the  rate  of  cooling  for  noriMl  alloy  solidification  causes 
the  rate  of  approach  to  equilibrium  solidification  to  be  slow;  too  slow 
for  the  extension  of  solidification  tiaies  to  be,  in  itself,  a  practical 
means  of  eliminating  pile roheterogenei ties  in  cast  structures. 

Copqiarison  of  the  solidification  curve  derived  from  the  therpwl 
data  of  the  slow  cooled  ingots  to  the  solidification  curve  calculated 
for  the  case  of  limited  solid  diffusion  indicated  that  the  latter  case 
does  apply  to  alloy  solidification.  In  addition,  piicrosegregation  measure¬ 
ments  indicated  that  there  was  an  increase  of  solute  within  the  solid 
phase  during  dendritic  growth  and  the  results  showed  good  agreement  with 
the  computations  for  the  conditions  of  limited  solid  diffusion. 
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I.  IHTROOUCTION 


This  report  suRMarlzes  c^ontlnulng  research  on  solidification  of 
alMinum  alloys  of  the  Casting  and  Solidification  Section  at  Massachusetts 
Institute  of  Technology  under  Army  Sponsorship  (now  Army  Materials  Research 
Agency)  through  Frankford  Arsenal.  In  recent  years  work  has  concentrated 
on  solidification  structure,  its  control,  its  relation  to  casting 
characteristics,  and  its  influence  on  mechanical  properties. 

Aieong  other  results  of  fundamental  or  practical  value,  techniques  have 
been  evolved  from  the  research,  which  are  now  employed  commercially,  to 
produce  quality"  castings;  i.e.,  castings  with  strengths  up  to 

twice  as  high  as  those  previously  available,  and  with  ductilities  up  to 
seven  times  greater. ^  These  improvements  in  properties  have  been  achieved 
in  existing  alloys,  primarily  by  control  of  (1)  microsegregation,  and 
(2)  microporosity. 

It  has  also  been  evident  that  more  complete  control  over  microsegre¬ 
gation  and  microporosity  (particularly  the  fornmr)  should  permit  development 
of  now,  wholly  different  and  substantially  stronger,  casting  alloys.  As  one 
example,  castings  with  tensile  strengths  of  70,000  psl  have  been  cast  in  our 
own  laboratory;  this  compares  with  50,000  psi  as  the  maximum  tensile 
strength  in  the  current  quality"  casting  specif ication.  ^  A 

research  program  is  now  underway  in  at  least  one  laboratory  to  develop 
alloys  which  possess  these  high  tensile  properties  and  have  adequate 

3 

ductility  for  engineering  applications. 


2. 


Oev«lopai«nt  of  wrought  aluminum  alloys  has  also  bean  limited  by 
microsegregation  and  microporosity  in  the  original  cast  ingots.  For 
example^  only  those  alloys  in  which  microsegregation  can  be  held  to  a 
required  minimum  can  be  readily  worked,  or  develop  adequate  mechanical 
properties  in  the  wrought  state.  An  example  of  the  as  yet  untapped 
potential  of  wrought  aluminum  is  that  a  recent  study  has  shown  that  with 
sufficiently  careful  processing  controls,  tensile  strengths  in  excess  of 

4 

100,000  psi  can  be  obtained  in  wrought  products  produced  from  cast  ingots. 

To  achieve  mechanical  properties  in  cast  or  wrought  products 
significantly  higher  than  those  now  available,  a  deeper  fundamental 
understanding  must  be  obtained  of  factors  affecting  such  structural 
variables  as  microsegregation  and  microporosity.  The  time  is  particularly 
ripe  for  such  study  since  (1)  a  variety  of  tools  can  now  be  employed  to 
aid  such  a  study  that  were  not  available  only  a  decade  ago,  including  the 
electron  microprobe  and  high  speed  computers,  and  (2)  a  decade  of 
experience  has  delineated  important  variables  and  shown  the  great 
improvement  in  properties  of  castings  and  wrought  material  that  is 
possible. 

This  report  summarizes  a  first  year's  effort  in  obtaining  fundamental 
information  as  discussed  above;  it  is  concerned  primarily  with  microsegre¬ 
gation.  Detailed  computer  analyses  of  the  solidification  process  are 
reported  and  comparisons  are  made  with  experiment.  Preliminary  suggestions 
are  given  as  to  how  the  principles  developed  may  be  applied  to  produce  more 
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castings  or  Ingots.  Work  following  that  raportad  hara  is 
l>alng  diractad  along  tha  dual  paths  of  Improving  insight  into  the  solidi¬ 
fication  procass  and  of  applying  Inforwation  davalopad^  on  a  laboratory 
scala,  to  altar  and  improva  cast  structures. 
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II.  ANALYSES  AND  COMPUTATIONS 

Once  the  trensport  and  equilibria  that  occur  on  the  microscopic  scale 
are  specified  or  assuined,  certain  macroscopic  properties  of  a  solidifying 
system  can  be  cor  outed . In  this  section  the  fraction  solid,  phase 
composition,  heat  content,  and  associated  derivatives  are  computed  for 
three  limiting  conditions:  (1)  equiliorium,  (2)  normal  non-equilibrium, 
and  (3)  limited  solid  diffusion.  The  aluminum-rich  portion  of  the 
aluminum-copper  alloy  system,  a  simple  eutectic  between  the  K,  0,  and 
liquid  phases  (Figure  I),  will  be  used  as  an  example  for  computation. 

The  analytical  expressions  derived,  however,  are  general  and  can  be 
applied  to  any  binary  alloy  system  whose  phase  diagram  is  accurately  known. 

Previously,  laborious  hand  calculations  were  necessary  to  carry  out 
this  type  of  calculation  for  any  but  alloy  systems  whose  phase  diagrams 
were  composed  of  straight  lines.  The  availability  of  an  IBM  7090  digital 
computer*  for  the  evaluation  of  the  expressions  made  feasible  carrying 
through  the  calculations  for  any  binary  alloy  system.  The  organization  of 
the  general  programs  for  computing  and  recording  the  results  of  the 
analyses  is  included  as  Appendix  B,  and  description  of  the  expressions 
fol lows . 

A.  Fraction  Solid 

For  a  solidifying  alloy  of  know  phase  diagram,  the  fraction  solid  at 
any  temperature  can  be  c::1culated  by  the  appropriate  application  of  mass 

*  This  work  was  carried  out  in  part  at  the  Computation  Center,  Cambridge 
Massachusetts . 
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b«laf>c«s.  An  alJoy  systeas  naving  «  ccwiponcnts  and  solidifying  through  a 
ragion  of  n  solid  phases  and  p  liquid  phasas  will  have  a-1  Independent 
sass  balances  of  the  fora 


-  (1) 


Cq  “  overall  weight  fraction  of  coea^nent  i. 


the  average  weight  fraction  of  consponent  1  In  the 
j  phase  of  the  solid  and  the  k  phase  of  the  liquid, 
respectively. 


the  weight  fraction  of  the  j  phase  of  the  solid  and 
k  phase  of  the  liquid,  respectively 


Consider  now  a  solidifying  alloy  of  a  specific  binary  system  (such  as 
that  of  figure  1).  At  any  point  in  the  two  phase  region  the  aass  balance 
for  solute  aiay  be  written  in  the  form 


vdiere  the  subscripts  i, 
phase,  one  liquid  phase. 


''0 


(2) 


j,  k  heve  been  deleted,  there  being  only  one  solid 
end  one  Independent  equation. 


Oif ferentlation  of  equation  (2)  gives  the  effect  of  en  Infinitesimal 
variation  in  the  system 

dfe^f.)  *  ^  ....  13) 


♦ 


The  symbols  used  in  this  report  are  tabulated  in  Appervdix  A 


o 


wnlcn  is  expanded  to 

As  iOTg  *5  t^'-e  SySter-a  is  closed 


■S  -  ^  ■  ’ 


(6) 


and  ec^atio’^  i^i  may  now  be  rewritte'"- 


(7) 


The  equality  exp'-essed  by  equation  (?)  is  depicted  for  an  inf initestaa^ 

a.moun:  of  sol  id!  f  icstlon,  dfj,  by  Figure  2a.  The  condition  before  the 

transformation  is  reqresented  by  the  1  * ne  “T  ...  and  the  condition  after 

c  o 

transformation  by  the  line  j  C ross “t atched  area  5  (represented  by  the 
left  side  of  equation  7)  must  equal  area  A  plus  area  C  (represented  by 
the  right  Side  of  equation  7) 


Equilibrium  Sol idlfication. 


Equilibrium  solidification  refers  to  a  liquid-solid  transformation 
that  occurs  with  complete  transport  within  both  liquid  and  sjlid  phases 
and  with  comoosltions  at  the  interface  strictly  as  predlctew  by  the  phase 
diagram.  The  compositions  of  the  solid  d  liquid  will  be  uniformly  those 
given  by  the  intersection  of  the  isothermal  at  that  temperature  (tie  line) 
with  the  solidus  and  liquidus,  respectively.  Figure  2b  depicts  the  change 


j(r 
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produced  by  an  anount  of  solidification  df^,  ^d  e<;uation  7  aodifi#>d  for 
this  special  cas«  is: 


. 

Cj*  «  waight  fraction  of  liquid  and  solid  prasant  at  intarfaca 
\*)  Straight  Phasa  Boundaries 


Ona 

solidus  and 


•athod  of  solution  jf  aquation  (8)  is  to  assuisc  that  both  the 
liquidus  are  straight  lines  described  by  and  where: 


a. 


where 


I 

!n 


i’ 


melting  point  of  solvent 

tao^ratura  of  invariant  transformation 

concentration  of  solvent  in  solid  and  liquid  phases  at  invariant. 


The  ratio  of  the  equilibrium  solid  concentration  to  the  equilibrium  liquid 
concentration  within  the  two  phase  region  is  a  constant  termed  the  equilibrium 
"partition  ratio". 


dCg*  •  kdCj^* 


01) 


and 


(12) 


3 


Mow  substituting  expressions  {11)  end  (12)  into  equation  (8),  separating 
variables,  and  L^oosing  1  im :  >  of  integration  yields 


1  1  -  (k-l)f. 


1 

(1  -  k) 


(13) 


Integration  leads  either  to  the  familiar  lever  rule 


f 


S 


(14) 


or  the  expression  used  for  coraputations 


S  (T  -  T)  (1  -  k) 
n 


(15) 


where  -  liquidus  temperature  for  alloy  Cq 
(b I  Curved  Phase  Boundaries 

Often  the  assunjpf'on  of  a  straight  liquidus  and  solidus  leads  to 
inaccurate  results  and/or  loss  of  information.  Such  would  be  true,  for 
example,  for  the  aluminum-zinc  system  (Figure  3)  or  the  aluminum-tin  system 
(Figure  4) , 

If  the  phase  boundaries  are  not  straight,  equation  (15)  no  longer  holds. 
Equation  (14),  however,  is  still  valid  and  may  be  applied  at  every  point  in 
.♦■he  path  Desj’ts  this  comp..,  tat  Ion  for  a  s8i  ics  or  a  i  uminum-copper 
alloys  of  differing  composition  are  plotted  in  Figure  5. 
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Although  th*  results  obtained  from  curved  line  calculation  often 
differ  appreciably  frc»  the  results  obtained  from  a  straight  line  approxi 
aiation  (equation  15),  equilibrium  solidification  is  rarely  approached  and 
the  significance  of  the  differences  will  be  discussed  in  relation  to  the 
non-equilibrium  calculations  (page  II  and  Figure  7) 


2.  Normal  Non-Equilibrium  Solidification. 


Normal  non-equi 1 ibriura  solidification  refers  to  the  condition  of 
(1)  complete  mass  transport  within  the  liquid  phase  (e.g.,  diffusion 
coefficient  is  infinite,  =  oo),  (2)  no  mass  transport  within  the  solid 
phase  {O5  *  0),  and  (3)  equilibrium  maintained  at  the  interface.  The 
solidification  assumptions  are  depicted  by  Figure  2c.  There  being  no  mass 
transport  within  the  solid  phase,  equation  (7)  is  simplified  by  the  dele¬ 
tion  of  the  terra  fs-df^  and  the  mass  balance  is  written  in  the  form 

-  Cs*)  dfj  -  (I  -  fj)  dCL* . (16) 

(a)  Straight  Phase  Boundaries 


Assuming  the  liquldus  and  solidus  to  be  straight  lines,  substituting 
exp'^essions  (ll)  and  (12)  into  equation  (16)  and  integrating  leads  to  the 
resul t 

/  Co 


f. 


I  - 


l/(l-k. 


.  .  (17) 


ut  liio  cApiession  used  for  computation 

/t. 


l/(l-k) 


1  - 


m 


.  (18) 


'm 
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(b)  Curved  Phas«  Boundaries 

When  the  phase  boundaries  are  not  sufficiently  straight  over  their 
entire  length  ot  justify  the  use  of  equations  (17)  and  (’.8),  the  phase 
boundaries  are  treated  as  a  series  of  straight  line  segments  that  follow 
the  actual  phase  boundary  within  the  experimental  accuracy.  Then,  for  a 
temperature,  T,  within  the  interval,  Tj  -  Tj.|,  that  the  phase  boundaries 
are  considered  straight,  define 


Cs 


J“1 


J-1 


(19) 


m, 


T,. 


J-1 


'Lj-1 


'-J 


-  .  .  (20) 


where  and  0$  are  weight  fraction  of  solute  at  intersections  of 
isothermal  at  Tj  (tie  line)  with  the  liquidus  and  solidus,  respectively, 
J  and  Csj_j  represent  similar  quantities  at  Tj_| .  Expressions  (19) 
and  (20)  combine  to  yield 


Cl*  -  C3* 


m, 


'"Sj 


c,‘  -  Cs.  ^ 


m, 


m. 


(21) 


Let 


"’L 


ms 


(22) 


cs 


(23) 


and  raplace  (C|_*  -  C^*)  tn  equation  (16)  by  AjC|_*  +  Bj 


dC^* 

AjC'l*  +  8j 


.  (24) 


Integrating* equation  (24)  ovsr  a  short  range  T*  to  T  within  the  interval 
Tj-i  to  Tj  using  as  lower  limits  of  integration  that  a  fraction  solid  f 5 ' 
exists  at  temperature  T'  and  liquid  composition  C|_'. 


f 


S 


(I  -  fj' 


)  ( 


Vl'  * 


+  B 


(25) 


To  calculate  the  value  of  f^  at  any  temperature  T  within  the  interval  Tj_j 
to  Tj,  equation  (25)  must  be  successively  evaluated  within  each  of  the 
intervals  from  the  liquidus  to  the  temperature  Tj_j  and  then  in  the 
interval  from  Tj_]  to  T. 

The  results  of  this  calculation  for  a  series  of  aluminum-copper  alloys 
are  plotted  in  Figure  6. 

In  Figure  7  results  obtained  for  several  systems  by  using  the  expression 
for  curved  phase  boundaries  (equation  25)  are  compared  to  results  obtained 
using  the  straight  line  assumption  (equation  18).  The  aluminum-iron  system, 
Figure  7b,  has  been  included  as  an  example  because  the  phase  diagram  is 
composed  oF  sfalght  lines  within  the  exper ir!w>nio?  accuracy  of  ine  determination, 
hence  results  calculated  from  both  expressions  coincide 

*  When  the  solidus  and  liquidus  are  parallel,  Aj  =  0.  Integration  then  gives 
instead  of  equation  (25),  f^  =  I  -  (1  -  fj')  exp  (Ci  '-Cl*/Bj) .  A  similar 
exception  exists  in  other  integrationj  fri  curved  phase  boundaries  that 
fol low. 


The  main  deviation  for  the  aluminum-copper  system.  Figure  7a,  is  that 
the  curved  boundary  calculation  predicts  a  higher  llquidus  and  thus  the 
curve  is  displaced  to  higher  temperatures.  The  predicted  amounts  of  non¬ 
equilibrium  second  phase  are  only  slightly  different. 

The  straight  line  approximation  is  not  even  qua  1 1 tative 1 y  applicable 
to  the  aluminum-tin  system,  Figure  7c. 

The  straight  line  approximation  is  valid  for  aluminum-zinc  alloys  in 
which  zinc  content  is  low.  Figure  7d,  in  spite  of  irregularities  in  the 
phase  diagram,  Figure  3-  This  is  due  to  the  fact  the  irregularities  occur 
at  higher  zinc  contents  and  the  low  zinc  alloys  are  mostly  solid  before 
entering  those  regions. 

3.  Limited  Diffusion  in  the  Solid. 

Limited  diffusion  in  the  solid  refers  to  the  assumptions  (I)  con^lete 
mass  transport  in  the  liquid,  (2)  mass  transport  in  the  solid  solely  by 
volume  diffusion,  and  (3)  equilibrium  maintained  at  the  interface  These 
solidification  conditions  are  depicted  in  Figure  2d  For  this  calculation 
an  additional  assumption  specifying  the  growth  morphology  is  taken  The 
growth  form  is  considered  to  be  plates  separated  by  a  distance  d.  The 
plates  start  to  grow  at  the  liquldus  temperature  and  complete  their  growth 
when  solidification  is  complete  at  a  time  later  (See  Flgu'-e  3  for 
growth  model  ) 
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An  increase  tn  the  solute  content  due  to  solid  diffusion  must  be  the 
result  of  a  flux  that  enters  the  solid  at  the  liquid-solid  interface.  Thus 
the  solute  increase  behind  the  interface  equal  to  the  flux  at 

the  interface. 

"  ■  '^i  “  "d" 

flux  at  interface  in  solid  phase  (gr/cm^) 
density  of  alloy  (gr/cm^) 
diffusion  coefficient  (cm^/sec) 

concentration  gradient  at  interface  {weight  fraction/cm) 
time  from  initiation  of  solidification  (seconds) 


where  » 

T  ■ 

0  = 

l‘i)  . 

dx/i 

8  = 


To  simplify  equation  (26)  the  amount  of  diffusion  is  assumed  relatively  small, 
i.e. , 


dCj* 

dx/ 


(27) 


(a)  Parabolic  Growth 


First,  let  the  rate  of  advance  of  the  Interface  be  parabolic  in  time 


X 


i 


=  X  V  8 


(28) 


where  Xj^  =  position  of  interface 


«  growth  rate  constant 


'k 


and 


d8  2Xi 

—  I  «  I  ■■■  a  -■  - . 

dXi  (>v)2  {X)2 


(23) 


For  th«  case  of  limited  solid  diffusion  and  plate-like  growth  form,  equation 
(7)  bacomes 


(Cl»  -  Cs*)  dfs  =  +  (1  -  fs)  dCL»  .  (30) 

where  ^  .  ...  (31) 


and  equation  (30)  may  only  be  applied  for  <"  I  due  to  the  assumption  (27). 

Straight  Phase  Boundaries.  Constant  Diffusion  Coefficient:  Once  again 
the  partition  ratio  is  assumed  to  be  constant  and,  in  addition,  the  diffusion 
coefficient  is  assumed  constant  over  the  solidification  interval.  Separating 
variables,  now 

dfj  I  dC^* 

kf;  -  fs  +  I  ■  Tl  -  it) 


and  integrating 


fs 


1  - 


(1  -  ^  k)/  (1-k) 


o<  k 


(33) 


or 


1  T  -  T 
/  _m _ L 

l^m  -  . 


(1-  »<  k)/(l-k) 


(3^+) 
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Curved  Phase  Boundaries,  Constant  Diffusion  Coefficient:  Equation  (30) 
can  be  evaluated  for  a  system  with  curved  phase  boundaries  by  dividing  the 
liquidus  and  solidus  into  a  series  of  straight  line  segments  (as  done 
previously  for  the  normal  non-equilibrium  case).  Using  definitions  (21), 
(22),  and  (23)  and  the  relation 


dCj*  -  (1  -  Aj)  dC^^ 


(35) 


equation  (30)  can  be  rewritten 


(AjC^*  +  8j)  dfs  -  (1  -  Aj)o<f^dCL*  +  (1  -  fg)  dCj_^ 


(36) 


Separating  variables  and  integrating 


''/Lj.  -  / 


where 


T-  (I  -  J  -  Aj]^  )  .  .  .  , 


( 


(37) 


(38) 


and 

f. 


1 

-  M  -  / 


S  I 


/Vl-  ^  ^  \  I  I 


Vl*  *  ^  ^ 


(39) 


Straight  Phase  Boundaries,  Temperature-Dependent  Diffusion  Coefficient 


The  diffusion  coefficient  of  solute  may  vary  by  more  than  an  order  of 
magnitude  over  the  temperature  range  of  solidification  of  an  alloy.  The 
alumlnum-i+  5  per  cept  copper  aHjy,  for  e?;amp’6,  ha^  a  solute  airrusion 
coefficient  of  I  62  x  10“9  cm^/sec  at  the  liquidus  and  at  the  eutectic  of 
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1.99  X  10"*®.  A  closer  approximation  is  made  by  allowing  the  diffusion 
coefficient  to  be  an  exponentially  varying  function  of  temperature  of  the 
form 

0  -  Oq  exp  (-Q/RTj.)  . (40) 


where  Oq  « 

Q  a 

R  - 

Now  cX 


empirical  constant 

heat  of  activation,  empirical 

gas  content 

absolute  temperature 

becomes  a  function  of  temperature  of  the  form 
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0<(T)  . 

(X) 


exp  {-(i/RTj^) 


(41) 


and  equation  (32)  is  integrated  over  short  temperature  ranges  between  the 
liquidus  temperature  and  the  temperature  of  consideration  At  each  interval 
the  proper  value  ofo<(  (T)  Is  substituted 


T“-'^(T)k 


[,  -  (I 


T  -  TI  ''  l-^''(T)k/  1-k  I 


Curved  Phase  Boundaries,  Temperature-Dependent  Diffusion  Coefficients: 

For  this  computation  the  dependence  of  0^  on  temperature  is  Inserted  into  the 
analysis  by  relating  o{  to  the  liquid  composition.  The  absolute  tempera¬ 
ture  at  some  position  within  the  interval  between  Tj_|  and  Tj  is  given  by 

\  (Clj  -  Cl>)  .  273  0 


(“*3) 


•nd  o<  (T)  Is  d*fln«d  «s  in  equation  (40)  Now  setting 


r (T)  =  I  -  (I  -  Aj)  (T)  (44) 


Use  for  each  successive  interetion: 


f  (T)fs*) 


/  A  C  ‘  +  B 

I  -i-k - 


i 

J 


Trans fortnet ion  at  the  Eutectic  Tereperature:  The  analysis  of  solidifi¬ 
cation  for  the  limited  solid  diffusion  case  has  carried  the  alloy  to  the 
point  of  eutectic  transforujation.  This  case  differs  from  the  previous 
two  cases  in  that  the  phase  redistribution  resulting  from  the  eutectic 
transformation  and  the  cooling  to  room  temperature  must  be  taken  into 
account.  In  the  equilibrium  esse  the  phase  distribution  is  simply  oDtained 
from  the  phase  diagram.  An  a  1 urainum-4 . 5  per  cent  copper  alloy  (Figure  1) 
is  totally  prlma-y  solid  soli-,  on,  K,  after  solidification  and  the  8  phase 
only  appears  as  the  aPoy  crosses  below  the  solvus  line.  In  the  normal 
non-equi i ibrium  case  there  can  be  no  mass  transport  within  the  solid  phase 
and  thus  no  diffusion  controlled  solid  state  change  can  occur.  An  alumlnum- 
4.5  per  cent  copper  alloy  that  nas  jO  '^rinary  sc<iid  ana  3  1%  liquid 
phase  just  above  the  eutectic  temperature  will  form  just  9.1^  eutectic 
This  requirement  is  relieved  for  tne  case  of  limited  solid  diffusion  An 
analysis  of  the  formation  of  the  primary  phase  at  t.ie  eutectic  temperature 
is  presented  below  An  analysis  of  solid  diffusion  below  the  eutectic  tempera¬ 
ture  will  be  presented  in  Section  0  (page  38) 
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Considering  the  eutectic  to  be  e  pseudophese  end  evaluating  the 
materials  balance,  equation  (1),  for  three  phase  equilibrium  of  a  binary 
alloy 

Cjf^  +  C^fg  +  C^fj^  “  Cq  {^(>) 

where  C^,  fj  r^fsr  to  the  prinaary  solid 
C^,  f^  refer  to  the  eutectic  solid 

eutectic  liquid 

for  an  infinitesimal  change 

d(Csfs)  +  +  d{Ci_\)  =  0  .  .  .  (^7) 

Expanding  equation  (47)  but  taking  into  account 

dCg  =  0;  dC^^  «  0,  °  h  " 

(48) 

df^  =  -df^  -  dfj_  .  (49) 

it  simplifies  to 

(Cs*  -  Cg*)  dfj  +  fsdCs  =  0  (49) 

The  second  term  of  equation  (50)  Is  evaluated  by  referring  to  equation 


fsdCs 


2Ps 

d 


(51) 


(26) 


Again;  tha  assuosptlon  is  takan  that  tha  solid  diffusion  dots  not  change 
significantly  the  concentration  profile  at  the  interface;  equation  (27); 


arvd 


fgdC^ 


40$  (Tg)  C^*  (i-k) 


d2  1  -  (!-<Kk)fs*(T^) 


:s2) 


where  refers  to  the  solid  fraction  just  above  the  eutectic 

teo^erature.  Rowrittlng  equation  (50)  and  setting  limits 


X  <S (Tf) 

/  / 


y 


4D3(T£)Cs*(I  -  k) 

d'2"[l  -  (1-  k)fsrE)j 

Og-f ’5(TE)2(d^Ax2) 


de 


.  (53) 


and  integrating 


fs(TE)  -  fs(Te) 


'  -  fslTj) 


I  -  (1  -  ^ 


(54) 


for  curved  phase  boundaries  the  ane'ogous  expression  is 
fs(Tc)  -  fs(7e)  '  -  f.'sdf/. 


(55) 


Evaluation  for  an  Aluminuni-4. 5  Copper  Alloy:  The  analysis  for  limited 
solid  diffusion  was  applied  to  an  aluminum-4.5  par  cent  copper  alloy  using 
the  published  phase  diagram  (Figure  1),  the  following  cc''Si;ants  in  the 
expression  for  the  diffusion  coefficient 
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0,  o  8.k  X  10"^  aTi^5«c*^ 

o 

Q  «  32,600  calories  mole”* 

and  a  parabolic  growth  rata  constant 

/  «  6  X  10“5  cm  sac”*^^ 

Tha  fraction  solid  at  each  temparatura,  l.a.,  tha  solidification  curva, 
for  tha  limited  solid  diffusion  co  dition  is  compared  to  those  for  tha 
equilibrium  and  norr<>al  -  nonaqui  1  ibrlum  limiting  cases  in  Figure  9*.  The 
calculated  curve,  based  on  equation  (45),  falls  between  the  lir.iting  curves 
for  tha  equilibrium  and  normal  non-equil Ibrlum  case,  the  position  of  tha 
curva  shifts  in  tha  direction  from  non-equilibrium  to  equilibrium  as  the 
value  of  X  decreases. 

For  parabolic  growth  of  dendrltlcal 1 y  freezing  alloys  7^  Is  equal  to 
one  half  the  dendrite  spacing  divided  by  the  square  root  of  solidification 
time;  i.e.,  from  equation  (28)  and  Figure  8 ; 

X  -  1/2  d//^  -  1/2^'  (36) 

where  0^  »  total  solidification  time 

-  empirical  ratio  of  dendrite  spacing  to  square  root  of 
solidification  time 

In  dendritic  solidification  of  ingot  cast  from  most  alloys,  \  is 

V* 

found  to  be  nearly  constant  over  a  considerable  range  of  freezing  From 


* 


Note’  Figure  9  Includes  some  experimental  data  described  in  Section  IV 
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a ) umlnum-coppsr  alloys  tha  valuas  of  ^  (*>2  A)  observad  ara  such  that  some 

solid  diffusion  Is  to  ba  axpactad  during  normal  solidification. 

(b)  Linaar  Growth 

Mow  consider  the  advance  of  tha  interface  to  be  linear  in  tin»c 


and 


-  u  0  . (57) 

dX. 

Ta  “  ^  . (58) 


For  tha  case  of  limited  solid  diffusion,  plate-like  dendrites,  and  linear 
growth,  equation  (7)  becomes 


-  Cs*)  dfs  •  ^  I  8Cs*  (’  -  fs)  ■  .  .  (59) 

where  "  20/ ud 


Tha  result  of  integration  for  tha  conditions  of  straight  line  phase 
boundaries  and  constant  diffusion  coefficient  is 


f 


S 


(1  +  cxf  j_k)  (I  - 


(61) 


Tha  analysis  has  also  been  applied  to  tha  conditions  of  curved  phase 
boundaries  and/or  tampcrature  dependent  diffusion  coefficient  in  a 
fashion  analogous  to  that  described  above  for  parabolic  growth. 


To  evaluate  the  result  at  the  eutectic  temperature  of  diffusion  in  the 
primary  solid  the  expression  for  f5d(;5  (equation  51)  is  evaluated  for  linear 
growth  rate 
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^Osd^)  (1  -  k)  Cs*  (T^ 
1  +o-<f  L  ^ 


Upon  substitution  for  this  term  in  the  materials  balance  for  the  eutectic 
temperature  (equation  50)  and  integrating 


f's(T£) 


■*  X  l'' 


UTLK-f'sfTE) 


(65) 


Similarly  for  curved  phase  boundaries: 


f  (Tj)  -  fjdj)  W,(l  -  Aj) 


I  -  f‘s(TE) 

1  *  .  Aj)  -  f5(Tj) 

...  (66) 


The  quantitative  results  of  this  analysis  for  linear  growth  rate  and 
limited  solid  diffusion  are  presented  in  the  next  section. 

If  the  growth  form  is  again  considered  to  be  dendrites  grown  from  the 
melt,  the  growth  rate  constant  u  is  half  dendrite  spacing  divided  by  the 
solidification  time;  l.e.,  from  equation  (57)  and  Figure  8. 

1  d 

u  “  -r  H — 

2  9f 

and 

1  cl2 

“  2  0~ 

where  still  is,  as  previously  defined,  the  empirical  ratio  of  dendrite 

spacing  to  the  square  root  of  solidification  time,  "^he  relation  of  the 
solidification  curve,  and  thus  the  microsegregation,  to  the  behavior  of 
receives  closer  attention  in  the  next  section  (page  29) 
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8.  Prop«rtlf$  of  th»  Sol  Idlf  lotion  Curve 

1.  Comparison  of  Several  Aluminum  Systems. 

The  solidification  curve  (temperature  versus  fraction  solid)  is 
important  in  the  prediction  of  many  casting  properties  and  the  interpreta¬ 
tion  of  certain  solidification  stuctures.  As  examples,  the  solidification 
curve  aids  in  understanding  relative  tendency  for  alloys  to  hot  tear,^® 
relative  shrinkage  characterisitics  of  different  alloys,'^  and  it  may 
affect  dendrite  morphology  as  well. 

The  solidification  curves  for  a  series  of  compositions  in  each  of  nine 
binary  aluminum  alloy  sy’^tems  were  computed.  The  results  for  the  case  of 
normal  non-egul 1 Ibrium  are  presented  in  Figure  10.  The  ordinate  values  in 
these  plots  are  not  degrees  as  previously  but  rather  a  dimensionless  tempera 
ture  scale  defined  by 


where  Tg  will  be  referred  to  as  the  reduced  temperature. 

For  an  alloy  that  transformed  uniformly  over  its  sol idif ication  range,  the 
curve  would  be  a  straight  line  of  slope  minus  one  (-1) 

Note  the  degree  of  similarity  that  exists  frcm  system  to  system  in 
Figure  10.  One  difference  between  the  curves  is  the  amount  of  solid  that 
forms  in  an  initial  temperature  interval  (e  g.,  the  first  five  per  cent 
t«nperature  decrease)  Comparing  equivalent  alloy  contents  of  the  systems 


a  I  uminuni~zinc ,  a  1  jmi nurn-nangnes iuf^ ,  and  a  1  unin jm-c opper  indicates  that 
the  an>ount  of  solid  to  form  in  the  same  drop  of  reduced  temperature 
increases  in  the  order  stated.  These  differences  may  be  expected  to 
influence  not  only  solute  redistribution  but  perhaps  also  dendrite 
morphology  in  the  different  systems. 

2.  Rate  of  Solid  Formation. 

Through  their  description  of  the  solidification  curve,  the  first 
derivative  of  the  fraction  solid  with  temperature,  df^/dT,  and  with  con¬ 
centration,  df^/dCj,  are  valuable  quantities  for  (!)  analyzing 
solidification  phenomena  and  (2)  insertion  in  solidification  calculations 
(e.g  ,  the  calculation  of  rate  of  shinkage  of  heat  '•ealease  that  follows). 
These  derivatives,  obtainable  by  rearrangement  of  the  differentia!  mass 
baJances  (equations  8,  16,  30,  and  59),  are  summarized  in  the  tabulation 
that  fol 1 ows 

The  reciprical  of  the  slope  at  the  liquidus  of  the  solidification 
curve,  (df 5/dTp) I ,  are  plotted  for  several  solute  contents  of  the 
aluminum-magnesium,  aluminum-copper,  and  aluminum-zinc  systems  in  Figure 
11  The  temperature  rate  of  solid  formation  decreases  sharply  with 
fraction  solid  (see  inset  Fig're  11)  Note  that  the  rate  of  solid  forma¬ 
tion  of  the  aluminum-zinc  system  is  about  five  fold  higher  than  the  two 
former  systems 


Case _ First  Derivative _ Straight  Line  Assumption _ Curved  Line  Assumption 
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3.  Composition  of  the  Solid. 

{«)  Composition  Within  the  Primary  Dendrite 

The  *B>ount  of  segregation  (coring)  which  results  in  each  of  the 
three  idealized  solidification  cases  discussed  in  Section  IIA  is  of 
interest.  The  minlm-jm  segregation  will  be  experienced  for  the  equilibrium 
case;  the  maximum  segregation  for  the  normal  non-equilibrium  case,  and  an 
intermediate  amount  of  segregation  for  the  limited  so. id  diffusion  case. 

As  an  extension  of  the  pre.ious  analyses  the  composition  at  the  interface 
and  the  average  composition  of  the  primary  dendrite  are  calculated  as 
fol lows . 


Interface  Composition  as  Function  of  Temperature:  The  composition 
present  at  the  interface  must,  for  each  solidification  case,  be  just  that 
given  by  the  phase  boundary  at  that  temperature.  For  straight  phase 
boundaries  interface  equilibrium  is  established  between  the  following 
compositions  in  each  case. 


and  also 


k  Cj_* 


_k_ 


(85) 

(86) 


(87) 


(88) 
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For  curved  phase  boundaries  tha  Interface  co(npo5ltions  take  the  form: 


C 


(89) 


and 


(I  -  Aj)  -  Bj 


(90) 


dC^* 

1 

dT 

dCs* 

1  -  Aj 

dT 

.  (91) 


(92) 


Interface  Ccwiposltion  as  a  Function  of  Fraction  Solid:  It  is  of 
interest  to  express  Cj*  as  a  function  of  fraction  solid,  as  this  relates 
closely  to  the  concentration  profile  within  a  dendrite.  The  relations, 
using  straight  phase  boundary  assumptions,  are  as  follows: 


Equilibrium:  Cj*  *> 

D 

-  (1-k)  fj] 

dCj* 

(1  -  k) 

^  2 

df,  ■ 

Cs* 

,  .k-1 

Normal  Non-Equilibrium: 

dCs*  Cs*  (1  -  k) 

d?r  ■  tt-TsI-  . 

Limited  Solid  Diffusion:  Cj*  -  kC^  (1  -  |j  kj 
(parabolic  growth) 


. (93) 

. (94) 

. (95) 

. (96) 

(k-l)/(l-o(k) 

fs) 

. (97) 
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dCs*  Cg*  {)  -  k) 

dfj  "  1  -  (l-o<k)f5 


(98) 


Limited  Solid  Diffusion: 
(linear  growth) 


C 


1  +  k  - 

i^r  - 

i  +  c./  J_k 


(99) 


dCs*  Cg*  (1  -  k) 

dfg  I  +  ^Vj^k  -  fjj 


(100) 


Average  Composition  of  the  Solid:  For  equilibrium  the  composition  within 
the  solid  is  uniformly  the  same  as  that  at  the  interface;  that  is,  the  average 
composition  is  equal  to  the  Interface  composition.  For  the  normal  non- 
equilibrium  case,  the  solid  composition  does  not  change  after  the  interface 
passes.  For  the  limited  solid  diffusion  case  the  solute  content  within  the 
solid  changes  continuously  and  approaches  the  interface  composition  For  each 
of  th©  three  cases,  the  average  composition  of  the  solid  at  any  point  during 
the  solidification  is  given  by  the  relation: 


C 


S 


(lOi) 


and 


(1  -  fg)  dC^* 

~17~  Tt" 


(102) 


where  the  derivatives  relating  to  the  appropriate  conditions  are  substituted 
on  the  right  side  of  equation  (102)  The  average  composition  as  a  function 
of  the  fraction  solidified  is  plotted  for  a  case  of  limited  solid  diffusion 
and  the  limiting  cases  for  an  aluminum-h  3%  copper  alloy  in  Figure  12. 
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Krora  the  solidification  curve  for  the  limited  solid  diffusion  conditions 
it  Is  only  f>ossib)e  to  calculate  the  interface  cocnposltion  and  the  average 
coflsposi tion .  The  computation  of  the  actual  segregation  pattern  at  each 
interval  of  solidification  is  computed  by  a  numerical  analysis  technique 
described  in  Section  D  (page  35) 

(b)  Amount  of  Second  Phase 

The  solidification  curves  for  several  compositions  of  tha  systems 
a  1  uminuta- copper ,  al umlnum '-magnesium,  and  alurainum-zinc  were  computed  using 
the  limited  solid  diffusion  conditions  and  several  values  of  ^  (=»d'^/9^) 

for  both  paraboli-c  and  llnea'  growth.  The  fraction  eutectic  formed  during 

tb  A 

solidification  as  a  function  of  is  plotted  for  several  compositions  of 
aluminum- cop per  having  parabolic  growth,  Figure  13,  and  linear  growth, 

Figure  14. 

For  the  a  1  u.mlnum -copper  alloys  up  to  the  maximum  solid  solubility 
limit  at  5.65  per  cent  copper  a  change  of  ^  (=2/.)  from  5  x  lO'^  cm/sec*^^ 

to  1  0  X  10*3  i5  sufficient  to  shift  the  behavior  from  essentially  equili- 
bri-ui  to  normal  non-equilibrium,  These  results,  which  are  summarized  by 
Figures  13  and  14,  are  supported  by  numerical  analysis  calculations  described 
later  (section  O) . 

The  results,  as  plotted  in  Figures  13  and  14.  show  the  sensitivity 
of  the  segregation  behavio,  of  several  aluminum -copper  alloys  to  .‘he  /alue 
of  (the  ratio  of  the  growth  eien>o»c  size,  e  g.,  the  dendrite  spacing,  to 
the  square  root  of  the  '•ol  idification  time)  Figure  15,  which  delineates 
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the  value  of  to  produce  a  given  aniount  of  non-equl  I  ibriura  eutectic,  is 
derived  by  a  transformation  of  Figure  13  A  horizontal  line  drawn  on  Figure 
13,  as  shown  at  one  per  cent  retained  eutectic,  intersects  the  curve  for 
each  coesposition  at  the  value  of  ^  just  necessary  to  form  a  given  jasount 

of  eutectic,  e.g.,  one  per  cent.  The  intersections  of  several  horizontal 
lines  are  plotted  in  Figure  15.  Thus,  the  area  below  the  curve  marked  0.5 
is  the  locus  of  points  for  which  one  half  per  cent  or  less  eutectic  will  be 
forced:  similarly,  for  the  curves  marked  2.5,  5.0,  10  0. 


The  results  for  the  aluminum  magnesium  alloys  up  to  the  maximum  solid 
solubility  limit  at  14.9  per  cent  magnesium  indicate  a  shift  in  V  from 

V 

1  X  10  ^  cfti/ssc^^^  to  2  X  10  ^  is  sufficient  to  shift  behavior  from  the 
non-egui  1  ibriura  type  to  essentially  equilibrium.  For  a  1  umtnura-zlnc  allovs 
up  to  10  per  cent  zinc  a  shift  from  ^  =  2  x  10”^  cm/sec’-'^^  to  ^  =>  2  x 

10"^  cm/sec*'^^  is  sufficient  to  cause  t  ge  in  solidification  bohavior. 

C ,  Heat  Content  of  a  Solidifying  Alloy 

The  detennination  of  the  solidification  curve  for  an  alloy  makes  possible 
calculation  of  a  variety  of  solidification  parameters,  the  actual  measurement 
of  which  is  often  extremely  difficult.  The  computation  of  two  such  groups  of 
porpe'‘tles  has  been  included  in  the  FORTRAN  program  fo*"  macroscopic  properties 
described  in  Appendix  B.  These  are- 


1.  The  specific  voluroe  of  the  phases  and  the  amount  and  --ate  of 


sh  r inkage . 
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2.  Th«  h«*t  content  of  the  phases  ard  the  rate  of  the  heat  release. 

Only  the  latter  coaputation  is  discussed  in  this  year's  report. 

The  total  heat  released  during  the  liquid  to  solid  transforeation  of 
a  binary  alloy  of  components  A  ard  S  depends  only  on  the  initial  and  final 
states  of  the  alloy.  The  heat  released  at  any  temperature  durin;  the 
transfoneation,  however,  is  intimately  dependent  on  the  solidification  path. 
This  dependence  on  the  path  is  significant  when  the  composition  changes  of 
the  solid  and  liquid  phases  during  solidification  are  considerable 

General:  The  heat  content  of  a  binary  alloy  at  any  temperature  in  a  two- 
phase  region  is  expressed  by  the  following  relations: 


Ht  *  ^$1^5  • 

where 

Hs  =  CjH'j  +  (1  -  C5)H"3  +  (T,  C.)  . 

Hl  -  +  (1  -  (T,  C^*) 

Having  used  the  following  expressions  for  the  solute 


(103) 

(104) 

(105) 

(106) 

(107) 
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Analogous  expressions  wore  u>ed  for  the  solvent  oarameters, 
where  single  primes  refer  to  solute  and  double  primes  refer  to  solvent  and 


Hy  =>  total  heat  content  of  alloy  relative  to  a  base 
temperature,  Tg  (calories/gram  alloy) 

Hs,  «  heat  content  of  solid  and  liquid  phases  (calories/ 
gram  phase) 


T 1  T” 


melting  point  of  pure  component  (’’C) 


C‘p,,  C'p  =  heat  capacity  of  pure  solute  in  the  solid  and  liquid 
^  states,  respectively  (ca 1 ories/gram°C ) 


.F« 


=  heat  of  fusion  of  pure  solute  (calories/gram) 


M  M 

H  (T,C^),  h[(T,C*)  -  heat  of  mixing  of  solid  and  liquid  phases  at 

temperature  T  and  average  compositions  and  C*, 
respectively  (calories/gram) 


The  assumption  has  been  made  that  the  properties  of  a  phase  are 
accurately  represented  by  its  average  composition.  This  is  nssentlaily 
use  of  the  Kopp-Neuroann  rule  of  heat  capacities. 


Equally  general  as  the  previous  expression;  are  the  following  derivates 
that  are  the  rate  of  heat  release  of  the  phases  and  the  aggregate. 


dHj  dHj  dM^  dfs 

df  “  dT~  dT“  ^  6T~ 


(108) 


where 


dHs  _  dCc 

—  -  CjC-p,  +  (l-Cs)C"p5  ^  (H'^.H-s)  57- 


dT 


-N  H  \  \  H  \ 

s  ,  ^  I 

\  C,  /  dT 
C  ^  ^  T 


(109) 
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dT 


C*C'  +  (1-C*,  )C''  +  (K',+H",  ) 


dC», 


Pl 


dT 


M  \ 

.k  , 


dC*. 


r;T 


C*. 


'*l/t 

(1 10) 


The  rate  of  he*t  release  with  increase  of  solid  fraction  (dHj/df^,  dHj_/df5, 
dH5/df5)  is  simply  obtained  by  multiplying  the  appropriate  derivative  above 
by  dT/dfj  (see  Table  l) . 


For  a  dand ri tica 1 1 y  freezing  alloy  the  liquid  is  considered  to  be 
homogeneous  on  the  microscale  and  thus  the  terms  that  depend  on  the 
solidification  path  are  f^,  dC^/dJ.  dT/dfg.  The  method  of  evaluation 
of  these  parameters  has  been  discussed  previously. 


The  analysts  has  been  applied  to  the  particular  case  of  an  aUsminum-i+  5 
per  cent  copper  alloy.  The  method  of  evaluating  each  of  the  thermal  ■'■ar9- 
mete-s  using  data  from  K.  K.  ^elley^^  and  from  0  Kubashaewskt  and 
J.  A.  Catterail'^  is  sunwarizud  in  Appendix  C. 

The  heat  content  is  plotted  as  a  function  of  temperature  Ir  Figure  16 
and  as  a  function  of  the  solid  fraction  in  Figure  17  The  hear  co''tent 
curves  of  Figure  16  are  shaped  similarly  to  the  so' idif ication  curves  for 
the  same  alloy  (Figure  9)  The  curves  for  the  non-equi • ib'lum  ease  differ, 
In  detail,  in  the  follc/wtng  manner, 

1.  The  curves  are  extended  in  the  direction  of  t’.«  abscissa  During 
sol  idif  icat.i.on  of  an  alloy,  the  sensible  heats  of  the  so’ id  and 
liquid  phases  tnust  be  removed,  In  addition  to  the  heat  of  fusion  of 
the  solidifying  fraction 


2.  The  ratio  of  heat  released  during  the  eutectic  hold  to  heat  released 
during  the  solidification  interval  is  not  equal  to  the  weight  frac¬ 
tion  eutectic  formed  During  the  eutectic  transformation,  which  is 
isothermal,  no  sensible  heat  is  removed  from  the  solid  or  liquid 
phases.  Only  neat  of  fusion  is  released  In  addition,  the  8  phase, 
which  is  one  of  the  eutectic  constituents,  has  an  appreciably  nega¬ 
tive  heat  of  formation,  approximately  twice  the  number  of  calories 
per  gram  as  the  negative  heat  of  nixing  of  the  eutectic  liquid 

3.  There  is  an  Inflection  in  the  heat  content  curve  near  70  per  cent 
solid.  The  composition  of  each  phase  and  of  the  transforming  fraction 
increases  continuously  during  solidification.  The  heat  capacity  of 
each  phase,  the  heat  of  mixing  of  the  liquid,  and  the  heat  of  fusion 
of  the  solidifying  fraction  arc  each  dependent  on  the  composition 
Reference  to  Figure  17  indicates  that  the  heat  that  need  be  removed 
for  each  fraction  solidified  increases  as  solidification  proceeds  and 
the  solute  content  increases  in  each  phase. 

If  the  heat  released  for  each  fraction  solidified  were  constant,  then  a 
heat  content  curve  could  be  transformed  to  a  solidification  curve  using  the 
construction  of  StonebrooK  and  Sicha.'^  Tne  slope  of  the  curve  in  Figure  17 
equals  the  amount  of  heat  that  need  be  extracted  to  form  a  gram  of  solid  at 
a  particular  stage  of  solidification.  It  is  clear  that  the  slope,  dHj/df^,  is 
not  a  constant  throughout  the  sol  idt<^ication  interval  of  an  aluminum-4.5  per 
cent  copper  alloy,  rather,  dHy/dfg,  continually  increases  as  solidification 

Hoveycr,  for  some  calculations  (e.g.,  qualitative  application  of 


progresses . 
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the  Stonebrooke  and  Sicha  calculation)  a  value  of  dHy/dfj  averaged  over  the 
solidification  interval  would  not  be  tto  crude  an  approximation 

0 .  Diffusion  Within  the  Dendrite 

I .  Ana  lysis. 

The  computation  of  the  limit  solid  diffusion  case  (Part  A“3,  page  12) 
delineates  the  conditions  necessary  for  significant  solid  diffusion  during 
sol tdif ication .  Enhancement  of  this  diffusion  reduces  the  severity  of 
microsegregation  and  is,  therefore  of  great  engineering  interest.  Previous 
calculations,  however,  do  not  fully  describe  microsegregation  in  that  they 
do  not  describe  the  distribution  of  solute  within  the  dendrite.  A  more 
complete  description  is  obtained  by  the  application  of  numerical  analysis 
techniques  to  the  diffusion  equation. 

The  model  used  in  the  numerical  analysis  computation  is  analogous  to 
the  model  of  Section  A-3  and  is  depicted  by  Figure  8.  The  pertinent  condi¬ 
tions  and  assumptions  are  described  below. 

1.  The  growth  forms  (or  dendrites)  are  considered  to  be  plates 

separated  at  their  centers  by  the  final  dendrite  spacing,  d.  The 
plates  start  to  grow  smoothly  at  the  liquidus  temperature  and  each 
plane  face  advances  a  distance  d/2  within  the  solidification  time 
0r. 

2  End  and  corner  e. fects  are  neglected 
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3  Th*  solut*  transport  within  the  solid  phase  is  solely  by  volume 
diffusion,  The  diffusion  coefficient  is  an  exponential  function 
of  temperature  described  by  equation  {^0}. 

4.  Density  differences  due  to  solidification  shrinkage  are  neglected. 

5.  Temperature  is  uniform  throughout  the  dendrite  plate  at  any  given 
time,  9  Although  a  casting  raay  be  solidifying  under  a  steep 
thermal  gradient,  the  temperature  difference  between  any  two 
dendrite  plates  will  always  be  small. 


6.  The  rate  of  advance  of  the  interface  is  taken  to  be  a  continuous 
function  of  time. 


parabolic  9  ....(111) 

or 

1  inear  «  u  9  ....(112) 


7.  (a)  In  the  first  analysis,  diffusion  behind  the  interface  is 

considered  not  to  alter  the  rate  of  growth.  For  this  assumption, 
the  interface  composition  at  each  position  of  the  moving  solid- 
liquid  interface  is  given  by  the  normal  non-equilibrium  expression 
for  interface  composition,  equation  (95),  which  becomes  for  the  co¬ 
ordinates  of  this  medel 
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s 


k-l 


KC^  (I  -  2X^/d) 


(113) 
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(b)  A  second,  mere  accurate  analysis,  but  one  tha^  requires 
more  ceraputer  time,  assumes  the  envelope  of  interface  cempesition 
is  given  by  the  limiKid  solid  diffusion  expressions  for  the  inter¬ 
face  composition,  ooi/4tien  (97)  for  the  parabolic  growth,  equation 
(99)  for  the  linear  growth.  Results  were  obtained  using  both 
analyses . 


Solid  Diffusion  Within  the  Sol idif icatien  Range:  Diffusion  within  the 
solid  phase  for  the  above  described  model  was  studi.^d  using  Pick's  second 
litw; 

jr  2/' 

(114) 


dd 


✓ 

_ 

5x2 


0 


with  the  initial  and  boundary  conditions  implied  by  the  model: 


when, 

T  «« 

then  0  ■  0,  X^  «  0,  and  »  kC^ 

.  .  (115) 
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X  - 

«  c*  . 

.  .  (116) 

(where  C*^  is  given  either  by  the  nortim!  non-equilibrium  or  the  limited 
solid  diffusion  expression  as  stated  in  assumptions  7  above). 

Using  the  method  of  finite  differences  equation  (114)  is  transformed  to 

EJ...k  -  .  i:j,k-l  -  Cj-l.k-ll 

^  A9  L  AX  AX  J 

.  .  .  .  (117) 

where  the  subscript  J  refers  to  steps  in  the  X  direction  ar.J  k  refers  to 
stops  In  time.  Solving  equation  (117)  for 
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(HB) 


where 


H 


DAS 


(119) 


Equatien  (118)  is  evaluated  Iteratively  over  J  and  k  and  the  boundary 
conditions  are  simulated  by  the  sequence  of  the  steps  taken  by  the 
computer  solution. 


Solid  Diffusion  During  the  Eutectic  Transformation:  Aqain,  as  the 
liquid  remaining  at  the  eutoctic  t^perature  is  transforming,  diffusion 
occurs  within  the  primary  solid.  This  diffusion  is  calculated  by  evaluating 
equation  (118)  at  all  the  J  positions  within  the  primary  solid  solution  for 
successive  iterations,  each  in  duration,  from  the  time  of  reaching  the 

outactic  temperature  te  time  9^  when  freezing  is  complete 

Diffusion  in  the  Solid  Solution  Region:  As  the  alloy  cools  frem  the 
eutectic  to  room  temperature,  a  certain  amount  of  homogenization  can  occur. 
If  this  homogenization  is  significant,  then  one  would  expect  to  see  less 
second  phase  and  more  level  concentration  gradients  at  room  temperature  than 
were  present  immediately  after  solidification  Once  the  time-temoerature 
path  for  the  alloy  is  determined,  the  above  analysis  can  be  extended  to 
calculating  the  amount  of  solid  state  homogenization.  Just  as  for  the 
solidification  range,  two  classes  of  cooling  curves  are  considered. 

The  condition  of  linear  advance  of  the  liquid  solid  interface,  equation 
(112),  for  the  model  requires  a  constant  rate  of  heat  removal,  q,  that  is 

q  =  constant  =  C|  ....  (120) 
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Such  <  heat  flew  conditien  wituld  be  largely  true  for  (1)  an  ingat  cast  sa 
that  interface  resistance  to  heat  flow  is  the  limiting  factor,  (2)  an  ingot 
cooled  by  the  "differential  thermocouple  technique"  described  later  in  this 
report,  and  (3)  for  "unidirectional ly"  solidified  ingots,  solidified  under 
steep  thermal  gradients,  for  distances  not  to  close  to  the  chill  and  tem¬ 
peratures  not  far  from  T£. 

If  the  heat  released  in  solidification  were  the  same  for  each  fraction 
solidified  and  equal  to  H‘  {cal ories/gram) ,  then 
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Nr  in  the  solid  region. 


^  d  dT 

q  “  "/Ss^  2  da 


(122) 


where  Cp^  =  heat  capacity  of  solid  (calories) 


Substituting  (122)  into  (121)  and  rearranging 
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Integrating 
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The  condition  of  parabolic  advance  of  the  liquid-solid  interface, 
equation  (111),  for  the  model  requires  a  parabolic  rat®  of  heat  removal. 


that  is. 
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(125) 


Such  a  heat  flow  conditian  would  be  largely  true  for  temperatures  nat 
to  far  below  Tg  when  resistance  ta  heat  flow  is  primarily  thermal  diffusivity 
of  the  maid  as  in  (1)  sand  castings,  (2)  chill  castings,  locations  near  the 
chill,  mold-metal  interface  resistance  negligible. 


For  parabolic  heat  removal. 
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and  again,  after  solidification 
(122)  into  (126),  rearranging. 


,  equation  (122)  will  hold, 
and  integrating 
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Substituting 


(127) 


Once  the  appropriate  time  temperature  path  has  been  selected,  equation 
(118)  can  be  evaluated  iteratively  from  the  eutectic  temperature  to  a 
temperature  at  which  the  diffusion  contribution  is  no  longer  significant. 


2.  Qualitative  Results;  Diffusion  During  Freezing. 


An  important  parameter  in  the  numerical  analysis  technique  is  M 
(defined  in  Equation  119).  It  is  a  characteristic  of  the  analysis  that 
when  X  represents  the  amount  of  solid  to  form  In  a  solidification  time 
.A  8 ;  for  H<<2  the  gradients  will  tend  to  levol  out  as  fast  as  the  dendrites 
grew;  for  M»2  there  will  be  little  diffusion  and  segregation  will  approach 
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the  taaxlnjum.  The  parameter  M  can  be  related  te  the  grewth  model  threugh'^ 
as  defined  first  in  equation  (56) 

d  -  . (128) 

It  is  deflwnstrated  belew  that,  within  the  assumptions  of  the  abeve 
analysis,  ever  any  range  that  the  value  ef  the  preportienality  factor 
for  an  alley  is  truly  a  constant,  the  solute  segregation  pattern  will  be 
the  sanoo  fo^  every  dendrite  spacing  and  coaling  rate.  The  enly  difference 
will  be  that  of  scale  (the  ratio  of  the  dendrite  spacings).  For  two  ingots 
n»ade  from  an  alloy  with  a  constant  value  of  ^  ,  the  solute  content  at 
dendrite  centers  will  be  identical,  the  solute  contents  at  equal  fractions 
of  the  dendrite  spacing  will  bo  identical,  and  the  amounts  of  non¬ 
equilibrium  second  phase  will  be  identical. 

Let  the  positien  of  the  interface  be  proportional  to  the  n^^  power  of 
the  solidification  time  (similar  to  equations  111  and  112), 


e"  . (129) 

A 

Then,  the  following  relation  exists  between  ^  and  pj 
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Consider  then  the  diffusion  in  the  slice  during  the  time  the  interface 
moves  from  Xj  =  (J  -  I )  da/  te  X2  =  J  d/a,  and 
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wh«re  a  »  number  of  spacas  taken  across  the  dendrite.  Then, 


0  =  00-8, 
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Now  evaluatina  H 


(134) 


where  *  2*'^'"'  -  (J  -  1)“^'  ^  is  a  numerical  factor 

dependent  on  the  relative  position  of  the  calculation  and  not  on  the 
absolute  distance.  Thus  M  is  dependent  only  on  D  and  ^  I'ut  not 
individually  on  d  or  In  particular  for  parabolic  growth,  n  =  1/2: 
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for  linear  growt.h,  n  =»  !  : 


M  = 


2a  D 


a 


(136) 


It  should  be  pointed  o>.  *  here  thjr  the  calculation  of  the  solidification 
curve  for  the  li-nited  solid  diffusion  case  also  gave  the  result  that  the 
amount  of  non-equi  i  tbriu.m  second  phase  was  depenOenr  only  ©n  2^  ^  *nd  not 
on  d  or  0^:  indivioua  •  1  y  Usirg  tiie  '•^lations 


Linear  growth  : 

Parabolic  growth-  ^  ^ 
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Equations  (34)  »ind  (61)  can  be  rewritten 


which  are  seen  to  depend  only  on  ^  ^ 


(139) 


.  .  (140) 


If  «  value  of  ^  is  known  for  an  alloy,  then  it  is  possible  to  casimate 
the  effectiveness  of  solid  state  diffusion  by  computing  20/^^. 

2 

If  •  »  20/^  ^  4  1  tendency  for  maxiroum  segregation 

2 

c><  ’  =  20/^  ^  I  significant  reduction  of  segregation 

2 

=  20/^  >  >.  1  tendency  for  minimurn  segregation 

Previous  researchers  have  reported  to  be  a  constant  over  a 

8  9 

considerable  range  of  dendritic  solidification  conditions.  '  Closer  study, 
however,  may  denxjnstrate  that  ^  is  not  a  constant  but  actually  a  function 
of  the  solidification  time.  Representing  the  relation  as 

. OM) 


when  m  “  0,  0  q  is  a  constant  and,  as  shown  above,  the  segregation 

pattern  will  not  alter  with  solidification  time. 


When  m^O,  the  value  ef  decreases  as  the  sel idif tcatien  titne 
increases  and  the  severity  ef  --^ring  will  be  lessened  fer  slew  ceellng 
rates  when  m^O,  the  value  ef  ,  Increase‘s  a?  the  sel idif icatien  time 
increases  and  cering  becernes  less  severe  fer  a  high  ceeling  rate. 

3.  Quantitative  Results;  Oiffusien  During  Freezing. 

The  numerical  analysis  technique  described  abeve  was  applied  te  an 
aluminum-4.5  per  cent  cepper  alley  assuming  a  parabelic  rate  ef  interface 
advance  (cendltien  111)  using  a  limited  solid  diffusien  envelepe  (assump¬ 
tion  7b)  and  using  the  fel lowing  numerical  value: 

0  - 1 

0^  -  0.084  cm  sec 

e 

0.  =  32,600  calories  mele"^ 

^  ^  «  1.2  X  10  ®  cm^sec  ^ 

The  concentration  profiles  at  progressive  stages  of  sel idif ication 
and  cooling  te  room  twnperature,  plotted  in  Figure  18,  demonstrate  the 
build-up  of  the  solute  in  the  solid  behind  the  advancing  interface.  Just 
abeve  the  eutectic  temperature,  the  solute  content  at  the  center  ef  the 
dendrite  is  predicted  at  1,32  per  cent  copper.  This  value  is  to  be  com¬ 
pared  te  the  0  61  per  cent  cepper  predicted  by  the  normal  nen-cqul 1 Ibr ium 
cenditions.  Just  below  the  eutectic  temperature,  the  selutc  content  at 
the  dendrite  center  is  predicted  at  1.34  per  cent  cepp'ir  and,  at  room 
temperature,  1.38  per  cent  copper. 
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The  fractien  s®lid  is  calculated  from  the  numerical  analysis  results 
by  subtracting  the  salute  centcr.t  as  integrated  ever  the  primary  phase 
frera  the  everall  alley  selutc  content  and  censidering  the  difference  to 
be  dispersed  in  the  secend  phase,  i.c.,  at  the  liquid  or  eutectic  coenposi- 
tien;  e.g.,  at  the  eutectic  temperature; 


For  the  results  plotted  in  Figure  18,  the  fractien  of  primary  solid 
selutien  just  above  the  eutectic  temperature  is  0.927;  just  below  the 
eutectic  temperature,  0.930;  at  reom  temperature,  0.935.  For  normal  non- 
oquil ibriura,  0.909  is  the  predicted  fraction  solid. 

The  amount  of  solute  build-up  within  the  primary  solid,  is  of  course, 

A 

sensitive  to  the  value  of  dependence  is  demonstrated  by  the  plot 

of  centerline  composition  versus  ,  of  Figure  19.  The  transition  from 

normal  non-equilibrium  behavior  to  equilibrium  behavior  occurs  in  the 
range  of  ^  ^  from  10"^  to  10”^*^  cm^sec"'.  The  range  of  the  transition 
predicted  by  the  numerical  analysis  technique  is  in  accord  with  the  range 
previously  calculated  by  the  analytic  expression  (45)  for  limited  solid 
diffusion  that  is  depicted  in  Figure  13. 

The  weight  per  cent  of  eutectic  predicted  at  the  eutectic  temperature, 
after  eutectic  transformation,  and  at  room  temperature  are  plotted  in 
Figure  20  as  a  function  of  )  ^  for  the  assumptions  (1)  parabolic  growth, 


he 


(2)  straight  line  phase  beundarles^  and  (3)  narmal  non~agui  I  Ibr iunt  envelape. 
It  is  clearly  demonstrated  by  the  curves  that  the  most  significant  reduction 
of  segregation  due  to  diffusion  in  the  primary  solid  occurs  during  solidifi¬ 
cation  and  not  during  direct  cooling  of  the  solid  casting. 

The  dotted  line  of  Figure  20  is  a  comparable  curve  calculated  by  the 
analytic  expression  for  limited  solid  diffusion,  parabolic  growth,  and 
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straight  line  phase  boundaries  (equation  42).  As  the  value  of  X  decreases, 
and  consequently  the  value  of  approaches  1,  the  analytic  expression  tends 
to  overestimate  the  amount  of  hfwnogenizatlon.  This  result  is  due  to  the 
assumption  stated  by  equation  (27)  that  solid  diffusion  will  net  change 
significantly  the  gradient  of  solute  at  the  interface.  The  change  that 
dees  occur  is  the  consequence  of  (1)  the  increased  rate  of  solid  formation 
and  (2)  the  buildup  of  the  solute  behind  the  Interface.  The  analytic 
expression  is  nonetheless  sufficiently  accurate  for  the  common  range  of 
values  of  ^  . 

Two  envelopes  were  used  to  express  the  interface  compositions,  the  narmal 
non-equilibrium  and  the  limited  solid  diffusion,  soe  assumption  7,  page  36. 
The  points  plotted  in  Figure  20  are  calculated  for  a  limited  solid  diffusion 
envelope  whereas  the  solid  curves  are  for  a  normal  non-equilibrium  envelope. 
The  former  envelope  is  expected  to  be  more  accurate  and  the  latter  more 
convenient;  however,  it  is  evident  from  Figuro  20  that  closely  similar 
results  are  obtained  using  the  two  envelopes. 
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III.  EXPERIMENTAL  APPARATUS  AND  PROCErjRES 


A.  General 

In  the  study  of  ordinary  ingots  and  sand  castings  it  is  difficult,  if 
not  impossible,  to  separate  the  interdependent  variables  that  effect  soli¬ 
dification.  It  is  expedient,  therefore,  to  utilize  equipment  and  procedures 
designed  to  independently  control  the  various  solidification  conditions  and 
isolate  the  important  parameters. 

Apparatus  used  for  the  controlled  solidification  studies  a^d  described 
below  (Figures  21  -  25)  was  designed  to  permit:  (1)  extended  solidification 
with  a  constant  rate  of  heat  removal,  (2)  programmed  time  temperature  paths 
during  solidification  and  subsequent  cooling,  (3)  unidirectional  cooling  with 
heat  removal  by  air  or  water,  (4)  mechanical  stirring  du''ing  cooling,  and 
(5)  mechanical  vibration  during  cooling. 

Specific  requirements  of  the  apparatus  were  that  it  would. 

(1)  Provide  a  wide  variety  of  solidification  times,  ranging  from  one  half 
hour  to  one  thousand  hours,  (2)  provide  an  environment  for  the  metal  that 
would  prevent  contamination  even  after  exceedingly  long  periods  of  contact, 

(3)  permit  close  control  and  observation  of  the  progress  of  solidification, 

(A)  incorporate  a  mold  of  such  a  shape  as  to  facilitate  mathematical  analysis, 
(5)  employ  a  relatively  insulating  mold  material  of  low  total  heat  content, 
and  (6)  permit  solidification  studies  to  be  made  with  a  minimum  set-up  time 
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1 .  Mel  ting  Practice. 

Charge  materials  for  al  uminun.-copper  heats  were  aluminum  (99.99+ 
purity)  and  electrolytic  copper.  The  metal  was  prepared  for  casting  in 
three  stages.  First,  copper  was  alloyed  with  liquid  aluminum  which  was 
held  and  stirred  above  1700°F  with  a  high  frequency  induction  furnace.  The 
resulting  hardener  alloy  of  nominal  composition  50  per  cent  aluminum-50  per 
cent  copper  was  cast  into  thin  ingots  in  cast  iron  pig  molds  Ne<t,  one 
hundred  sixty  pound  heats  of  aIuminum-4.5  per  cent  copper  sl  o'’'  wee  made 
by  combining  the  virgin  aluminum  and  the  50-50  hardener  alloy  in  a  gas-fired 
air  furnace.  The  resulting  alloy  was  cast  into  approxi-iiately  6  pound  slab 
ingots  in  cast  iron  melds.  Finally,  individual  heats  of  twenty  to  forty 
pounds  were  made  by  remeltin,  the  master  ■»’loy,  stirring  f#r  sixty  seconds, 
degassing  with  chlorine  for  sixty  seconds,  and  then  pouring  into  the  test 
ingets.  Clay-graphite  crucibles  were  used  for  all  beats.  Ail  meta'  tools 
and  pig  .molds  were  coated  with  either  fiberfrax  or  zircon  wash  and  dried 
before  use  The  procedure  described  above  was  used  to  assure  complete 
dissolution  and  uniform  distribution  of  the  copper,  consistent  composition 
from  ingot  t©  ingot,  and  minimal  contamination  from  iron  a-'d  silicon 

2.  The  Mold. 

A  cylindrical  steel  mold  was  fabricated  by  we'ding  a  */8-inch  steel 
plate  to  the  bottom  of  a  seamless  steel  tube  of  l/8-incb  wa) '  thickness. 

Two  mold  sizes  were  employed,  5"inches  diameter  and  2-i,''2  inc'^es  diameter 
respectively,  each  approximately  10  inches  in  height  Mo'id-metal  reaction 
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was  praventad  by  caatlng  the  inside  maid  surface  with  Fibarfrax  cement* 
and/ar  Fiberfrax  paper.  The  top  ef  the  maid  was  cavered  by  a  transite 
plate  The  insulating  ability  af  the  maid  was  adjusted  by  wrapping  layers 
af  Fibarfrax  paper  araund  the  outside  surface  of  the  steel  shell. 

3.  Solidification  Furnace,  Centrals,  Thermal  Measurements. 

The  furance  far  the  controlled  solidification  studies  was  a 
resistance-wound  recirculating  air  furnace  capable  of  a  maximum  temperaturo 
of  ISSO^F,  Figures  21  -  25.  Two  Brown  proportional  band  control  units  were 
available;  ene  a  duration  adjusting  type  (DAT),  the  second  a  current 
adjusting  type  (CAT).  The  calibration  accuracy  of  the  controllers  is  +0.2 
per  cent  ef  full  scale  (1500‘’F  far  OAT,  300°F  far  CAT).  Two  Speedamax  N 
model  S  recorders  having  interchangeable  ranges  and  scales  sf  one  millivolt, 
ten  millivolts,  and  fifty  millivolts  were  also  available  far  thermocouple 
reading  and  recording.  .A  constant  millivoltage  was  connected  in  opposition 
to  the  thermocouple  when  the  scale  range  in  use  was  less  than  the  thermo¬ 
couple  output. 

Thermocouples  wore  fabricated  from  Hoskins  special  chromel -a lumel  wire 
with  a  guaranteed  accuracy  of  +3/8  per  cent.  The  wire  was  ca'ibratod  before 
shipment  and  the  calibration  was  checked  In  the  laboratory  in  the  range  200'  F 
te  1220°F  If  the  rate  of  temperature  change  was  2  F/mi'-  »r  'esj,  ^8  gauge 
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wtr«  was  used  and  th«  thtrmacaupl®  baads  caatad  with  fiberfrax  The  caating 
liraited  Interactian  with  the  melt  and  raducad  stray  currant  pickup.  If  a 
faster  raspansa  ware  raqulrad,  uncaated  22  gauga  tharmacauple  wira  was  used. 

B.  Procedure 


t 


) 


» 
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In  the  fallawlng  paragraphs,  dataliad  praceduras  are  described  that  were 
amplayed  far  praducing  tast  castings  in  the  central  led  sal idif icatian  appara¬ 
tus  in  the  fallawing  ways!  (1)  extended  sal idif icatian,  (2)  interrupted 
sal idif icatian,  (3)  isatharmal  sal idif icatian,  (4)  pragrafwnad  sei idif icatian, 

(5)  unldlrectianal  sal idif icatian,  (6)  sal idif icatian  with  stirring,  and 
(?)  sal idif icatian  with  vibratien. 

In  this  repart,  hawaver,  detailed  results  are  presented  anly  far  *-nose 
experiments  which  fall  under  1  and  2  abava  (plus  same  limited  data  f -am  5) 
Experi.ra'iis  in  the  remaining  categaries  were  largely  explaratary  in  nature 
during  ^he  fiscal  year  cavared  by  this  rapart,  these  experiments  are  cantinuing 
and  will  be  repartad  in  detail  in  the  next  annual  repart  af  this  canfact 

1  Extended  So! idif Lcation 

Ingats  were  selidified  aver  extended  periods  ranging  from  eight  hours 
to  ana  thousand  hours  while  maintaining  (I)  the  rata  af  heat  extractian 
constant  and  (2)  the  temperature  gradient  in  the  melt  very  shal'aw  (lass  than 
0.1°F  inch).  The  experimental  setup  is  sketched  in  Figure  2!  and  the  proce¬ 


dure  is  summarized  below 
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1.  Ing«t  ro*1ds  ware  placed  In  the  recirculating  air  furnace  and 
preheated  at  1300®F  fer  at  least  twelve  heurs. 

2.  The  melds  were  filled,  and  the  melt  temperature  stabilized  at  the 
furnace  temperature. 

3.  The  melt  and  furnace  were  ceeled  with  centre!  responding  to  a 
differentia]  thermeceuple  that  measured  the  difference  between  the 
temperature  ef  the  melt  and  the  temperature  of  the  furnace  chamber. 

As  depicted  by  Figure  21,  the  differential  thermaceupl i  consisted  ef 
one  chremel -alumel  junction  in  the  melt  connected  in  parallel  with  a 
second  junction  in  the  furnace.  The  potential  across  the  differential 
thermeceupl a,  which  represented  the  temperature  difference  between  the 
two  measuring  points,  was  input  to  the  controller  which  acted  to  keep 
the  temperature  difference  at  a  sot,  constant  value. 

4.  When  the  recorded  melt  temperature  indicated  the  ingot  has  passed 
through  the  eutectic  hold,  the  ingot  was  removed  from  the  furnace  and 
water  quenched. 

The  rate  of  heat  extraction  is  directly  proportional  to  the  temperature 
difference  and  is  a  constant  as  long  as  the  difference  is  a  constant  The 
rate  of  heat  extraction  is  reduced  (and  the  solidification  time  extended) 
by  (a)  reducing  the  temperature  difference  or  (b)  increasing  the  insulation 
around  the  mold.  The  temperature  gradient  in  the  melt  is  simultaneously 
reduced  by  either  of  these  two  procedures. 
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At  tii»*s  as  many  as  six  individual  thermocoupl as  wara  used  in  a 
diffarantial  thermocouple  package.  Three  in  the  furnace  were  connected  in 
series,  thus  multiplying  the  millivoltage  reading  by  three.  Similarly, 
three  thermocouples  in  the  melt  were  connected  in  series.  When  these 
thermocouples  were  connected  to  form  a  differential  thermocouple,  the 
voltage  across  the  package  corresponded  to  three  time  the  actual  tempera¬ 
ture  difference.  Not  only  does  this  nf>ethod  scale  small  temperature 
differences  up  to  a  poin,  where  they  are  read  .nore  easily  and  contro’led 
more  accurately,  but  this  method  also  averages  temperatufe  variation?  in 
the  furnace  and  averages  error  in  thermocouple  readings 

2.  Interrupted  Solidification. 

A  variation  of  the  above  techniques  was  to  slow  cool  severe!  2-1/2 
inch  diameter  ingots  in  the  furnace  setting  the  differential  thermocouple 
control  from  just  one  of  the  ingots.  At  each  of  a  series  of  temperatures 
below  the  liquidus,  one  of  the  ingots  would  be  removed  from  the  furnace. 

The  mold  was  then  quenched  to  150°F  by  dropping  into  a  large  tank  of  water 
Less  than  two  minutes  elapsed  from  the  time  the  slow  solidification  was 
interrupted  and  the  metal  temperature  brought  to  150°F.  Thus  the  portion 
of  the  alloy  that  was  still  liquid  when  the  slow  solidification  was  ir^ter- 
rupted  was  cooled,  by  the  quench,  about  3600  times  faster  than  the  remainder 
of  the  al 1 oy . 


A  secord  technique  employed  to  interrupt  solidification  consisted  of 
pouring  lead  of  the  same  temperature  as  the  melt  into  the  mold  The  fraction 
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of  th«  alloy  still  liquid  would  ba  displaced,  but  the  solid  crystals  would 
be  trapped  by  a  fiberglass  screening  fastened  across  the  centnr  of  the 
laold . 


3.  Isothermal  Solidification. 

A  mold  with  a  minimum  amount  of  insulation  was  placed  in  the  furnace 
and  preheated  to  1300®F.  The  mold  was  filled  with  the  furnace  door  open  and 
the  metal  ten^erature  was  lowered  to  a  position  intermediate  lo  the  solidus 
and  the  eutectic.  The  solidification  temperature  w«s  maintained  For  a 
specified  time  and  then  the  mold  removed  from  the  furnace  and  quenched 

At  times  several  2-1/2  inch  diameter  ingots  were  cooled  to  the  solidifi¬ 
cation  twnperature  simul taneousl y  end  each  was  held  at  that  temperature  for 
a  different  length  of  time  before  quenching. 

I4.  Programmed  Solidification. 

A  Beck  Programmer  was  used  in  conjunction  with  the  control  units  to 
program  temperature-time  cycles  for  the  metal  and/or  furnace.  Figure  22  is 
a  sketch  of  one  experimental  set-up  which  was  employed.  The  programmer 
outputs  a  time-varying  millivoltage  that  either  added  or  opposed  t^e  output 
of  a  control  thermocouple.  Tne  algebraic  sum  of  the  two  mi  1 1 ivo 1 tages  was 
then  connected  to  a  control  unit,  which  maintained  it  at  a  constant  value 
Thus,  as  the  output  from  the  progranner  varied  the  contro.ler  changed  the 
power  input  to  the  furnace  to  the  point  where  the  control  thermocoup' es 
compensated  for  the  change.  The  programmer  is  capable  of  producing  up  to 
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thr««  diff»r«nt  linear  cooling  rates  within  the  specified  millivo  1  tage 
interval s . 

5.  Unidirectional  Solidification. 

A  chill  was  fabricated  connecting  four  1/4-inch  diameter  channels 
using  two  1/2-inch  channels  in  a  steel  block  5"  x  5"  x  1/2"  and  welding  it 
to  the  bottofls  of  a  5“inch  diameter  cylindrical  mold.  The  chill  was  coated 
with  Fiberfrax  cement  and  the  inside  and  the  outside  surface  of  the  mold 
walls  were  wrapped  with  Fiberfrax  paper.  The  experimental  set-up  is  sketched 
in  Figure  23-  Typically,  the  ensemble  was  placed  in  the  furnace  and  preheated 
to  1300°F.  After  the  mold  was  filled  and  the  metal  temperavjre  stabilized, 
either  water  or  compressed  air  was  run  through  the  chill  block  at  a  hi^;h  rate 
(while  the  furnace  temperature  was  maintained  above  the  liquidus).  Essen¬ 
tially  all  heat  flow  in  the  ingot  was  "unidirectional  1 y",  downward. 

6.  Solidification  with  Stirring. 

A  1/3  horsepowder  motor  and  1/2-inch  chuck  were  mounted  on  a  frame 
attached  to  the  top  of  the  furnace.  The  circuit  bringing  power  to  the  stirrer 
motor  contained  a  potentiometer  to  adjust  the  power  input  and  voltmeter  and 
ammeter  to  measure  the  power  input.  Various  paddle  arrangements  could  be 
inserted  into  the  chuck  and  one  arrangement  is  shown  in  Figure  24.  The 
paddle  and  shaft  were  generally  made  of  type  31&  stainless  steel  and  coated 


with  aluminum  and/or  fiberfrax. 
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7.  Solidification  with  Vibration. 

An  Eriex  Hi  VI  model  V3A30N  electro-magnetic  vibration  unit  having  a 
maximum  power  input  of  40  watts  was  mountad  on  a  frame  below  the  furnace 
and  coupled  to  the  mold  by  a  stainless  steel  rod  which  extended  through  the 
furnace  bottom  A  stainless  steel  plate  was  threaded  and  locked  onto  the 
rod  and  the  mold  itself  was  securely  bolted  to  the  plate.  Figure  25  is  a 
photograph  of  an  experimental  setup.  The  vibrator  power  could  be  turned 
on  at  any  time  during  solidification  and  the  vibrator  power  adjusted  by 


means  of  a  rheostat. 
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IV  RESULTS  OF  CONTROLLED  SOLIDIFICATION  EXPERIMENTS 

A  Thorwal  Analysis  of  Extended  Solidification  Heats 

Sev«ril  ingots  of  aluminum-4.5  per  cent  copper  were  solidified  over 
extended  periods  ranging  from  twelve  to  1000  hours,  at  a  constant  rate  of 
heat  extraction  using  the  differential  thermocouple  tecl.nique.  Photographs 
of  polished  and  etched  macrosections  of  ingots  solidified  twelve,  sixty,  one 
hundred  ind  one  thousand  hours  are  shown  in  Figure  26.  Thermal  record  of 
ingots  solidified  over  twelve,  thirty  six,  and  eighty  two  hours  are  plotted 
in  Figure  27.  The  plots  are  similar  in  shape  to  the  heat  content  (Figure 
16)  and  solidification  curves  (Figure  9)  calculated  earlier  for  this  alloy. 

The  slope  decreases  sharply  as  the  first  solid  appears  and  remains  relatively 
shallow  and  linear  for  a  good  part  of  the  solidification  time  The  slope  near 
the  end  of  solidification  is  close  to  the  value  previous  to  tHe  transformation . 
The  horizontal  portion  represents  the  eutectic  hold  and  final ‘v  the  slope 
increases  to  a  value  ch-iracteris tic  of  a  single  phase  field.  Neithe”  ti'e 
aajount  of  undercooling  below  the  liquidus  prior  to  the  initiation  o^  freezing 
nor  the  undercooling  below  the  eutectic  prior  to  the  initiation  of  the 
eutectic  transformation  was  eve.'  measured  as  greater  thar  two  degrees  Farepheit 

1  The  Solidification  Curve,  Comparison  with  Computations 

Due  to  the  fact  the  ingots  were  coo’ed  with  «-  constarr  'ate  of  heat 
extraction,  the  abscissa  of  the  thermal  record  differs  from  the  heat  zoi'tcnt 
by  a  constant  factor  and  the  coo'ing  curves  may  be  compa-'ed  :o  solidification 
curve  (fraction  solid  versus  temperature.  Figure  9),  using  the  method 
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Ston*brook  and  Sicha.'^  The  method  assumes  that  the  heat  of  fusion  liberated 
for  each  per  cent  of  transformation  is  the  same  and  that  the  heat  capacity  of 
the  solid  and  liquid  ar«  the  same.  The  heat  content  calculations  of  Section  II- 
C  demonstrate  that  these  assumptions  are  individually  inaccurate,  especially 
so  when  an  intermetal  1 ic,  such  as  CuAl2  is  one  of  the  transformation  products. 
But,  the  construction,  when  applied  to  an  alumlnum-4.5  per  cent  copper  alloy  is 
accurate  within  15  per  cant  during  the  bulk  of  solidification.  The  converted 
thermal  data  is  plotted  in  Figure  9  with  the  calculated  curves  for  the  three 
cases.  The  data  is  seen  to  agree  reasonably  well  with  the  case  of  limited 
solid  diffusion  and  a  value  of  X  ^  •  1.2  x  10"®. 


2.  The  Amount  of  Non-Equilibrium  Second  Phase. 


It  is  interesting  to  note  that  there  is  only  a  very  slow  increase  in 
the  amount  of  solid  diffusion  when  the  time  of  solidification  is  increased 
over  two  orders  of  magnitude.  The  ratio  of  the  time  of  the  eutectic  hold  to 
the  total  solidification  time  may  be,  through  the  calculations  of  Section  II- 
C,  equilibrated  with  the  weight  fraction  of  eutectic  formed.  These  values, 
summarized  in  Table  II,  change  only  slightly  in  the  runs  twelve  to  eighty  two 
hours.  There  was  no  hold  recorded  for  the  one  thousand  hour  heat,  however, 
some  eutectic  is  observed  in  the  nsacrostructure. 


It  is  obvious  from  the  macrostructures  (Figure  26)  that  the  dendrite 
spacing  is  increasing  along  with  the  solidification  time  Even  though  the 
time  for  solid  diffusion  to  take  place  may  be  increased,  the  dif<^usion  paths 
becomes  extended  and  thus  the  concentration  gradients,  rapresentati'.e  of  the 
driving  force  for  diffusion,  become  shallow  The  somewhat  simplified  mode! 
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TABLE  II 


Solidification 

Time  (hours) 

Duration  of  Eutectic 

Hold  (hours) 

Weight  Fraction 
Eutectic  (calculated*) 

12 

0.25 

.03 

36 

0.7 

.03 

60 

1.3 

.03 

82 

1.7 

.03 

1000 

None  Observed 

- 

•n4  c*lculations  of  Section  I  (in  particular,  equation  141)  indicate  that  if 
the  dendrite  spacing  increases  in  a  manner  such  that  the  ratio  of  the  square 
of  the  dendrite  spacing  to  the  solidification  time  (  remains  constant, 

the  segregation  pattern  wi!'  not  change.  The  experimental  evidence  clearly 
indicates  that  the  extension  of  solidification  times  is  not,  in  itself,  a 
practical  means  of  eliminating  microsegregation  in  aluminum-if  5  per  cent 
al loys . 


*  The  method  of  calculating  a  heat  content  curve  for  the  solidification 
interval  is  described  in  some  detail  in  Section  II-C :  heat  content  of 
solidifying  alloy.  The  fraction  eutectic  is  derived  from  the  calculated 
heat  content  curve  that  has  the  ratio  heat  released  at  eutectic'  total 
heat  released  within  solidification  interval  equal  to  the  ratio  from  the 
cooling  curve  of  time  at  eutectic  temperature:  time  within  solidification 
interval . 


59. 


B.  Hicros»qr«q«tion  H«sur«m€nts* 

1.  Hicroprobe  Measurements;  Unidirectional  Ingots. 

The  segregation  pattern  within  the  dendrites  of  a  series  of  ingots 
cast  from  aluminum-4.5  per  cent  copper  was  analyzed  quantitatively  with  the 
aid  of  an  electron  mi(:roprob6.  The  photomicrograph  of  Figure  28  is  a  section 
taken  perpendicular  to  the  heat  flow  direction  and  seven  inches  from  the 
chill  from  an  ingot  that  was  water  chilled  and  unidirectional ly  solidified. 
The  sample  was  prepared  for  analysis  by  first  polishing  meta 1 1 ographical 1 y 
and  then  placing  r  xcrohardness  indentations  to  mark  the  area  of  observation. 
The  electron  probe  traces  are  apparent  in  the  photomicrograph  and  the 
results  of  one  trace  are  plotted  on  the  same  figure  with  the  scale  of  the 
abscissa  such  that  horizontal  distances  are  the  same.  The  composition  along 
the  centerline  of  the  dendrite  element  Is  fairly  uniform  and  for  this  sample 
a  minimum  value  of  1.3  per  cent  copper  was  obtained. 

The  minimum  copper  content  measured  at  the  centerline  of  dendrite 
elements  of  samples  taken  from  locations  on  this  and  other  ingots  is 
summarized  in  Table  III  The  minimum  solute  content  expected  if  the  soli- 
diflcatien  followed  normal  non-equilibrium  cenditions  is  0  6  per  cent  copper 
for  an  aluminum-4. 5  per  cent  copper  alloy.  The  fact  that  the  minimum 
compositions  observed  are  about  1.3  per  cent  copper  indicates  that  the 


*  The  Ingets  used  for  these  microsegregation  measurements  were  prepared  by 
Or.  T.  F.  Bower  and  the  electron  microprobe  measurements  were  made  under 
his  direction. 
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TABLE  III 

A.  Unii<lr«ctiTial  I Y  Solidifl«^  Ingots 


Method  of  Chil 1 

Distance 

from  Chil ' 

Minimum  Copper 

1 . 

Water  Chi  1 1 

1/16 

inch 

1  .4 

3 

inches 

1 .4 

7 

inches 

1.3 

2. 

Air  Chi  1 1 

1/16 

inch 

2  2 

7 

Inches 

1.3 

3. 

Air  then  Water 

1/2 

inch 

0.8 

3/4 

inch 

0  9 

7/8 

1  h 

B. 

Slow  Cooled  Inqets 

1 . 

Interrupted  at 

120I°F 

5% 

Solid 

0  1 

2. 

Interrupted  at 

1  187‘'F 

40% 

So  '  id 

\  0 

•  ssumpttons  of  this  limiting  cas«  art  not  co(np>etei>y  valirf  and  rhat  on*  of 
th«  two  following  conditions  contributes  to  tho  discrepancy  (a)  thee  is 
a  significant  build-up  of  soiut*  of  the  inltia  gr*wth  forms 

due  to  limited  liquid  diffusion  or  (b)  some  so’utio-'.l  zing  occurs  Ir  the 
prlmaiy  dendrites  during  so  1  Idl  f  icar  ior  t^>*r*  a  bui'd-ip  of  so'jte 

in  the  liquid,  the  minimum  solute  co"t«nt  measu'ed  wi  tnln  tb*  soHd  phase 
would  be  the  same  throughout  so  i  Idl  f  ica  r  io”  On  ot'^er  hand,  solid 
diffusion  in  the  dendrites  would  cause  tre  mlnoKm  •oiL.te  to  g'^aduaily 
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incre«s«  during  s.«l  i^lf icatl»n  slmilcr  t»  th«  pr««iictl»n  of  Figure  18. 

Tw«  €xp«rifB«ntal  techniques  ef  interrupting  sel  i^if icatien,  described 
belc¥^,  were  used  te  deraenstrete  that  the  selute  centent  within  the  dendrites 
dees,  indeed,  increase  during  sel idificetien,  end  that,  therefere,  the 
increase  in  capper  cancentretian  at  the  center  af  dendrite  arms  abeve  the 
value  kC^  results  frere  diffusian  in  the  salid.* 

2  Hlcraprabe  Hcasu reman ts.  Interrupted  Sel id ificatian. 

First,  an  ingat  was  salidified  unidirectienal 1 y  by  air  chilling  until, 
accarding  ta  thermal  data,  the  mast  advanced  sal  id  raached  3/^  inch  Inta  the 
ingat.  Than  watar  was  put  thraugh  tha  chill  ta  effectively  quench  in  the 
existing  structure  The  dete  in  Table  I  shews  that  the  pri.-aary  set  id  phase 
hed  8  -  .9  per  cent  capper  at  the  time  af  quenching  The  regian  abave  'i/k 
Inch  had  a  nsrmal  minimura  selute  cantent  far  a  water  c.hilled  ingot,  1.^  per 
cent  cepper. 

Secend ,  ingats  were  slaw  caalad  ta  tamperatures  belaw  the  liquldus 
and  than  lead  af  tha  sees)#  tarnparature  es  the  alley  was  paured  into  the  ingat 
ta  displace  the  liquid.  In  this  case  <n  Inget  that  had  its  sal idif icatien 
interrupted  at  120l°F,  equivalent  te  5  per  cent  sel id,  hed  a  selute  centent 
af  0.7  per  cent  capper  within  tha  priinary  dandritss.  An  ingat  that  had  its 
sal idif icatian  interruptad  at  llS/^F,  equivalent  te  40  per  cant  s»!id,  had 
«  salute  content  af  1  0  per  cent  capper 

*  Alse,  relatively  simple  caiculatians  and  cemputer  calculation  now  underway 
as  a  part  ef  this  pragr«s  shew  that  under  cendltiens  af  sal idif Icatien 
censiderad  herein  there  can  be  na  significant  build-up  af  salute  in  the 
liquid  in  front  ef,  er  between  dendrite  arms. 
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Th.  .xp.rl«.„t.l  r.jults  subst.ntl.l.  th.  r.sults  .f  c.lcul.tl.ns 
which  inJlot.  th.t  llritei  S.IIJ  Wlffusl^,  can  .ccur  during  n.r».l 
s.lutflctun.  Th.  v.lu.  .f  th.  mlnlmu.  j.lut.  cnt.nts  »..sur..  Is 


c«n5lst.nt  with  .  v.lu.  .f  V  ^ 
•f  S«ctl«n  D  an^  Flgura  19. 


1.2  X  10”®  Crt 


fT  tha  calcu- 


i!, 
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V  CONCLUSIONS 

This  r«p*rt  covers  th«  first  y««rs  w#rlc  •n  a  cantinuing  prajact  far 
tha  funWanantal  stu^y  af  tha  sali^lf Icatlan  bahaviar  af  aluninuts  allays. 

An  analytic  stu^y  has  baan  af  micrasagragatian  in  Ingats  cast  fraos 

substantially  allaya^  matals  Tha  sal i^if icatian  pracass  has  baan 

iaaaliza^^  far  tha  saka  af  analysis,  by  tha  salectlan  af  the  mlcrascapic 
transpart  can^itians  an^  af  a  grawth  fDa<jel  with  tr<a  fallawing  results 

1.  The  sol Idif Ication  curve  (fraction  solid  as  a  function  of  temperature) 
for  binary  alloys  has  been  described  and  quantitative  results  ha% • 
been  presented  for  several  alutainuin  base  alloys  This  metnod  is 
easily  applied  to  ternary  (and  more  cotaplex)  alloy  systems  provided 
adequate  phase  diagram  information  is  available 

2.  Afjalyses  given  account  for  limited  diffusion  in  tne  solid  pnase  during 
solidification  and  during  cooling  to  room  temperature  Tnesc  analyses 
allow  more  accurate  description  of  the  solidification  conditions  and 
more  useful  analysis  of  the  solidification  structure  than  conp^ tat  ions 
based  on  the  usual  assumptions,  (a)  complte  diffusion  in  the  solid 
(equilibrium},  or  (b)  no  diffusion  in  tne  solid  (normal  non-equi I ibrium) 

3.  Tne  amount  of  non-equi 1 ibrium  phases  and  coring  present  in  the  final 

sol idif ication  structure  is  shown  to  be  dependen  on  tne  ratio  of  the 
square  of  the  size  of  the  growth  element  (for  example,  dend.^ite  arm 
spacing)  to  the  sol  idi  fication  time,  expressed  as  ^  For  the 

dendritic  solidification  of  an  aluminum-4  5  per  cent  copper  alloy  a 
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vaiwCC’"  nessc''ce.  ''eSoitinf'e  T.axI'Tjm 

segregation,  a  valjc  of  >.  ^  =  I  x  10”*^  >ou!d,  in  essence,  result  in 

an  eguiliDriuKj  structure.  A  value  of  ^  ranging  between  these 
limits  would  result  in  varying  amounts  of  solid  diffusion  and  c_rves  t^at 
describe  the  amount  of  non -«qu i 1 ibr turn  eutectic  formed  in  aluminum- 

\  O 

copper  alloys  as  a  function  of  '  have  been  oresanted 

U.  As  a  consequence  of  solid  diffusion,  the  solute  content  within  the 
primary  solid  will  increase  during  the  solidification  process  A 
numerical  analysis  technique  has  been  described  wnich  comoutes  the 
concentration  profile  at  various  solid  fractions.  The  analysis  applied 
to  an  aluminum-i*.  5  per  cent  copper  alloy  demonstrates: 

(a)  the  dependence  of  the  microsegregation  on  ■  ^ 

(b)  the  copper  concentration  of  the  primary  solid  phase  is  quite  uniform 
during  the  initial  stages  of  sol idi f ication 

(c)  whenever  there  is  substantial  solid  diffusion  t^'t  reduction  of 
microsegregation  during  cooling  in  the  solidification  range  will 
be  more  significant  then  the  leveling  tnat  will  occur  during  the 
eutectic  hold  and  during  cooling  of  the  soP  '  structure  to  room 
tenperatu  re 

5.  Once  the  ' o! idif ication  curve  of  an  alloy  has  been  computed  the  behavior 
within  the  solidification  range  of  many  physical  properties  can  be  des¬ 
cribed  witn  the  use  of  data  for  the  alloy  that  nas  been  measured  witnin 
the  single  phase  regions  As  an  example,  expressions  rave  been  derived 
for  the  heat  content  and  rate  of  heat  release  of  a  binary  alloy  .n  the 
sol idif ication  region  The  expressions  have  been  evaluated  for  an 
aluminum-U.5  Per  cent  copper  alloy 
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Through  the  enelysis  of  Ingots  cast  froti?  aluminum-4. 5  per  cent  copper 

alloy  axpertraental  verification  was  obtained  for  the  major  assumptions  of 

the  analyses. 

1.  Thermal  analysis  of  ingots  soliilfie^  by  a  differential  thermocouple 
technique  indicate  that  the  solidification  curve  of  aluminum-4.5  per 
cent  alloy  is  more  aptly  described  by  the  limited  solid  diffusion 
conditions  than  either  the  equilibrium  or  normal  non-equilibrium 
conditions . 

2.  Electron  microprobe  measurements  of  the  copper  distribution  within 
dendrites  from  unidirectional ly  solidified  ingots  showed  a  mininxim 
value  of  1.3  ■  1.4  per  cent  copper  rather  than  the  0.6  per  cent 
predicted  by  the  phase  diagram.  Examination  of  dendrites  from  ingots 
that  had  their  cooling  interrupted  at  various  stages  of  sel idif ication 
demonstrated  that  initial  solid  to  solidify  did  have  a  composition 
close  to  0.6  per  cent  and  that  the  solute  content  of  the  solid  increased 
during  the  solidification  process 

3  Observation  of  ingots  cooled  with  solidification  times  ranging  from 
twelve  hours  to  one  thousand  hours  demonstrated  that  the  dendrite 
structure  coarsens  with  increasing  solidification  tin>e  at  such  a  rate 
that  the  resulting  decrease  in  concentration  gradient  hinders  solid 
diffusion  to  an  extent  comparable  to  the  enhancement  brought  about  by 
the  Increased  time  aballable  for  diffusion.  The  reduction  of  micre- 
segregatien  resulting  from  extended  solidification  times  Is  too  slight 
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to  D«,  in  itself,  a  practical  means  for  eliminating  ne  nor-eqji ’ icrium 
micr®-inhoo>ogeneities  from  sol  i«  if  ication  structures 
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Figure  22;  Schematic  diagram  of  a  possible  experimental  set  up  for 
solidification  by  a  programmed  cooling  cycle. 
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unidirectional  sol  id  1 f 1  cal  1 un  of  an  Ingot. 
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Figure  25:  Photograph  of  apparatus  used  for  vibration  of  Ingots. 
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frequency  constant  for  diffusion 

volume  diffusion  coef'^icient  of  solute  in  solid 

weight  fraction  solid,  liquid,  and  eutectic, 
respectively 

heat  released  for  eacn  fraction  solidified 


heat  of  fusion 


heat  content  of  the  solid  and  liquid  phase, 
respectively 

heat  of  mixing  of  ln«  solid  and  liquid  phases, 
respectively 

heat  content  of  a  binary  alloy 

flux  at  the  interface  in  the  solid  phase 

partition  ratio 

finite  difference  melho-  •  s  ant, 


slopes  of  Uquldus  and  solidus 

slopes  of  liquidus  and  solidus  in  the  interval 

heat  of  activation  of  diffusion 

rate  of  heat  transfer 

gas  constant 

temperature 

base  temperature  for  heat  content  measurements 

temperature  of  an  invariant  reaction,  e  g  ,  eutectic 
t ransforiTiation  temperature 

temperature  in  degrees  keivcn 

liquidus  temperature 

reference  temperature  of  interval  of  phase  dicgram 

melting  point  of  pure  solvent 

growth  rate  constant,  linear  growth 

distance  along  growth  element 

position  of  solid-liquid  interface 

20/Uf^d 

empirical  ratio  of  square  of  dendrite  spacing  to 
solidification  time,  d^/8j. 

1  -  (1  -  Aj) 

growth  rate  constant,  parabolic  growth 
density  of  alloy 

tinve  since  beginning  of  solidification 


total  solidification  time 
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APPENDIX  B:  Organization  of  Computer  Programs 

Evaluation  of  the  solidification  parameters  Introduced  In  Section 
II  has  been  achieved  with  two  FORTRAN  programs.  The  first,  MACRO, 
coinputes  the  paraownters  that  are  descriptive  of  an  entire  phase  or 
the  entire  alloy  and  are,  at  least  hypothetically,  ammenable  to  macro¬ 
scopic  measurement  (pyroinetry,  calorimetry,  dilatometry,  etc.).  These 
are  the  parameters  discussed  In  Sections  II-A,  II-B,  II-C.  The  second 
program,  MICRO,  computes  the  concentration  profile  in  the  growth 
element  on  the  basis  of  microscopic  diffusion  conditions  as  described 
in  Section  II-D.  The  organization  of  each  is  described  by  a  flow 
diagram  of  the  executive  routine  of  each;  MACRO,  Figure  29,  MICRO, 

Figure  30, 

MACRO  consists  of  an  executive  sub-program  (MAIN)  that  calls  s 
sequence  of  computational  and  input/output  subroutines.  The  calling 
sequence  of  (MAIN)  is  presented  in  Figure  29.  Many  variations  have 
been  compiled  for  each  computational  and  output  subroutine  and  the 
particular  series  of  subroutines  selected  to  buildup  a  deck  for  a 
computer  run  is  a  choice  dependent  on  the  nature  of  the  alloy  system 
(e.g.,  real  or  Irrwigirary) ,  the  assumed  solidification  conditions,  and 
the  parameters  to  be  evaluated  and/or  compared.  The  general  purpose 
of  each  subroutine  is  indicated  in  Figure  29.  A  program  of  this  design 
was  deer?>ed  most  suitable  for  a  research  project  in  which  the  assumptions, 
parameters,  and  systesns  taken  under  consideration  are  continually 
changing  and  expanding. 


102 


MICRO  consits  of  an  executive  sub-prograw  (MAIN)  and  several 
computational  and  input-output  subroutines.  The  calling  sequence  of 
(main)  Is  delineated  by  the  flow  diagram,  Figure  30,  as  are  the 
functions  of  each  subroutine.  The  scope  of  job  to  be  run  on  the 
computer  is  read  in  as  part  of  the  data  and  thus  with  the  exception 
of  subroutine  CHICO  only  one  subroutine  of  each  name  is  needed 
Because  output  is  a  relatively  slow  computer  step,  the  output  routines 
have  been  designed  to  be  easily  altered  by  octal  correction  cards  when¬ 
ever  computer  time  must  be  conserved 
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APPENDIX  C:  Thermal  Data  for  Heat  Content  Analysis 
of  A1  -  4.5  Per  Cent  Cu 


I .  thermal  Properties  of  the  Pure  Elements 

K  K.  Kelley^^  fits  the  temperature  variation  of  heat  content  of  an 
element  to  an  expression  of  the  following  general  form 

H(T)  »  aT  +  bT^  +  C 

The  evaluation  of  this  expression  for  liquid  and  solid  aluminum  and 
copper  are 

-  |'(4.94)T  +  (1.43  X  !0-3)T^  +  (-l,604)j  7  (26.98)  (cal/gr) 

-  ^(5.41)T  +  (0.75  X  IO'3)t2  +  (-1,680)]  7  (63.54) 

Hl'  -  [(7.00)T  +  (330)]  7  (26.98) 

-  •[(7.50)T  +  (-20)]  7(63.54) 

and  the  heat  capacity  expression  is  the  fit  to  the  derivative  of  the 
above  expressions,  that  is 

Cp  -  a  +  2b  T 

and  Cp^’s  "  [(^-94)  +  (2.96  x  I0"3)T  |  7  (26.98)  (cal7gr‘’C) 

Cp^^'s  “  {(S.'+l)  +  (1.50  X  !0-3)Tj  7  (63.54) 

Cp^^L  "  (7.00)7  (26.98) 


(7.50)7  (63.54) 


2.  Heat  of  Mixing  Data 


(a)  Assume  to  be  negligible.  The  average  copper  content  of  the 
primary  solid  is  less  than  two  per  cent  (less  than  one  atomic  per  cent). 

(b)  From  0,  Kubashewski  and  J.  A.  Catterall^^  the  heat  of  mixing  of 
the  liquid  up  to  33  weight  per  cent  fits  the  relation 


hJ_  “  -{9000)Nq^  (cal /gram) 


where  Nr  *  atom  fraction  copper 


(26.98) 


Cl  (26.98)  +  0  -  C*l)  (63.54) 


Define  VD  such  that 


-  i  .  NcJ63.5‘>)  .  (1 


grams  liquid 

^Cu' (28.98)  =  gr*(n  atom  liquid 


Then 


(9000)  Nq^j  (cal/gr  liquid) 


To  evaluate  dH|^/dCL 


=  -(9000) 


+  N, 


^ 

"Cu 


which  becomes 


63.54 

9,000 


(c)  For  the  interiT>etan  tc  phase  CuAKCS) 


H  / 


-  3,250  (cal/gr  atom) 


and 


Hq  -  -  3,250  (cal/gr) 


tA 


